Symmetry protected fractional Chern insulators and fractional topological insulators 



o 

(N 

a* 
< 
o 

(N 



03 



i 

o 
o 



> 

(N 
O 

OV 

O! 



Yuan-Ming Lu 1 and Ying Ran 1 

'Department of Physics, Boston College, Chestnut Hill, Massachusetts 02467, USA 

(Dated: April 23, 2012) 

In this paper we construct fully symmetric wavefunctions for the spin-polarized fractional Chern 
insulators (FCI) and time-reversal-invariant fractional topological insulators (FTI) in two dimensions 
using the parton approach. We show that the lattice symmetry gives rise to many different FCI 
and FTI phases even with the same filling fraction v (and the same quantized Hall conductance a xy 
in FCI case). They have different symmetry-protected topological orders, which are characterized 
by different projective symmetry groups. We mainly focus on FCI phases which are realized in 
a partially filled band with Chern number one. The low-energy gauge groups of a generic a xy = 
1/m • e 2 /h FCI wavefunctions can be either SU(m) or the discrete group Z m , and in the latter 
case the associated low-energy physics are described by Chern- Simons-Higgs theories. We use our 
construction to compute the ground state degeneracy. Examples of FCI/FTI wavefunctions on 
honeycomb lattice and checkerboard lattice are explicitly given. Possible non-Abelian FCI phases 
which may be realized in a partially filled band with Chern number two are discussed. Generic 
FTI wavefunctions in the absence of spin conservation are also presented whose low-energy gauge 
groups can be either SU(m) x SU(m) or Z m x Z m . The constructed wavefunctions also set up the 
framework for future variational Monte Carlo simulations. 

PACS numbers: 71.27.+a, 71.70.Ej, 73.43.-f 



I. INTRODUCTION 

Phases of matters in condensed matter systems can al- 
most always be characterized by the Landau-Ginzburg 
symmetry breaking theory^. Experimental discovery of 
integer and fractional quantum Hall states in 2-D electron 
gas under a strong external magnetic field^ has provided 
striking counter examples of this paradigm. The frac- 
tional quantum Hall liquids arc particularly fascinating 
in the sense that their low energy excitations are quasi- 
particles carrying fractional electric charge^ and obey- 
ing anyonic statistics^. Although these liquid phases do 
not break physical symmetries, they are still different 
quantum phases. One measurable difference is their edge 
states: despite the fact that these liquids are all insula- 
tors in the bulk, they all possess certain edge metallic 
modest. In general different bulk phases host different 
edge states which can be detected by various experimen- 
tal probes such as electric transport^. 

A few years after the experimental discovery of inte- 
ger quantum Hall effcct(IQHE), Haldanc showed that the 
essence of it is not the external magnetic field^, by explic- 
itly writing down a lattice model Hamiltonian of IQHE 
with zero net magnetic field. However, it takes more 
than two decades for people to show that similar state- 
ment is true even for FQHE. Recently results from a se- 
ries of model studies^—, including convincing evidences 
from exact diagonalizations^— , indicate that fractional 
quantum hall states exist in the ground states of inter- 
acting lattice models, in the absence of an external mag- 
netic field. It is found that the ground state is likely 
to respect the full lattice symmetry. Here we call these 
fractional ground states spin-polarized "fractional Chern 
insulators" (FCI) to distinguish from the traditional frac- 
tional quantum Hall states in an external magnetic field. 



These proposed lattice models share a common feature: 
a partially filled nearly flat two-dimensional band with 
non-trivial band topology. 

The concept of band topology originates from the 
well-known TKNN index (or Chern number) of an IQH 
insulator—. In the past few years, this concept has been 
generalized to time-reversal symmetric systems, and trig- 
gers the theoretical and experimental discoveries of topo- 
logical insulators in spin-orbital coupled compounds in 
both two and three spatial dimensions^— . In two di- 
mension (2D), a time-reversal symmetric band insulator 
is characterized by a Zi topological index. Experimen- 
tally, HgTe quantum heterostructure has been shown to 
be a 2D topological insulator—. In the simplest limit, 
2D topological insulator can be viewed as a direct prod- 
uct of the up-spin and down-spin wavefunctions hosting 
opposite TKNN index. 

It is then quite natural to ask whether similar timc- 
reversal-invariant (TRI) versions of two-dimensional frac- 
tional topological insulators (FTI) exist or not, and there 
has been a lot of interest in this issued— . In the sim- 
plest limit when spin along the z-direction is conserved, it 
can be understood as the direct product of wavefunctions 
of the up spin and the down spin with opposite FQHE. 
Clearly this direct product is a fully gapped stable phase. 
In addition it must have non-trivial ground state degen- 
eracy on a torus even in the presence of a small S z conser- 
vation breaking perturbation, because the ground state 
degeneracy cannot be lifted by an arbitrary local per- 
turbation. So there is no question that in principle this 
fractionalized phase could exist. One important issue is 
whether this phase hosts stable gapless edge excitations. 
This problem has been studied in by Levin and Sterol. 
Another important open question is that whether TRI 
FTI can exist in a reasonable Hamiltonian. 
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In order to realize the FCI or FTI phases in experi- 
ments, one should find a compound with a nearly flat 
topological non-trivial band so that correlation effect is 
strong. Naively this is unnatural because usually a flat 
band is realized by spatially localized orbitals which do 
not support topological non-trivial hopping terms. How- 
ever, a very recent theoretical investigation^ on transi- 
tion metal oxide heterostructures indicates that a nearly 
flat topological non-trivial band can be naturally realized 
in the e g orbital double-layer perovskite grown along the 
[111] direction. Exact diagonalization in the same work 
shows that fractional quantum hall state can be realized 
in principle when the nearly flat band is partially filled. 
Because the temperature scale of the FCI/FTI physics in 
this system is controlled by short-range Coulomb inter- 
action, it can be a high-temperature effect. 

Fractional quantum Hall states, especially the non- 
Abelian ones, have been shown to be very useful build- 
ing blocks of quantum computers. If high tempera- 
ture FCI/FTI physics can be realized experimentally, 
it will certainly have deep impact in condensed mat- 
ter physics, including the efforts on topological quantum 
computation^. 

Motivated by the recent progresses on FCI/FTI 
physics, in this paper we try to address several impor- 
tant issues: what are the many-particle wavefunctions of 
2-D FCIs/FTIs? Can there be more than one FCI/FTI 
phases with the same filling fraction? If the answer is 
positive, can we classify these quantum phases ( or ground 
state wavefunctions)? 

Historically, Laughlin's wavefunctions of FQH states 
in a magnetic field 5 - have been shown to be one of the 
most important theoretical progresses in many-particle 
physics. It allows people to understand a lot of properties 
of FQH liquids in a compact fashion, including the frac- 
tionalized quasi-particle excitations^, topological ground 
state degeneracie s 30 ' 31 , as well as constructing the low en- 
ergy effective theories 3 ^—. Here in the case of FCI/FTI 
systems, analytical understanding of the ground state 
wavefunctions will help us extract various measurable in- 
formation in a similar way. 

Recently there was an interesting work to construct 
FCI wavefunctions by proposing a one-to-one mapping 
between the lattice problem and the magnetic field 
problem 3 ^. We would like to emphasize that the wave- 
function problem for FCI is related to that for the mag- 
netic field case, yet they are very different from each 
other. This is because the lattice symmetry of FCI is fun- 
damentally different from the continuum case of the 2-D 
electron gas. In fact, the recently discovered FCI states 
preserve all the lattice point group symmetry as well as 
translational symmetry^ 3 - Here in this paper, we point 
out that as a consequence of the lattice symmetry, there 
exist many different quantum FCI phases, all respecting 
the full lattice symmetry, even at the same filling fraction 
with the same quantum Hall conductance. These differ- 
ent FCI phases are distinct in the bulk in a more subtle 
way. One hand-waving statement is that the bulk quasi- 



particle excitations of these phases carry different lattice 
quantum numbers. These distinct FCI phases cannot be 
adiabatically connected with each other without a phase 
transition while the lattice symmetry is respected. Sim- 
ilar phenomena of distinct topologically ordered phases 
protected by symmetry is known in the context of quan- 
tum spin liquids 3 ^ and other low dimensional topological 
phases 3 -!. 

Now we outline the content of this paper. We start 
with the spin-polarized FCI at filling v = — (m is an odd 
number). In section [TT] the SU(m) parton construction 
of the fractional quantum Hall states (or spin-polarized 
FCI states) is introduced on a lattice, which is a natural 
generalization of the continuum cas o 5 ' 38 . We argue that a 
general FCI wavefunction could break the SU(m) gauge 
group down to Z m , and consequently the low-energy 
dynamics is described by Chern-Simons-Higgs theories. 
We explicitly write down the form of the electronic FCI 
wavefunctions which will be useful for future variational 
Monte Carlo study. We construct quasiparticle excita- 
tions of such FCI states. To demonstrate how lattice 
symmetry restricts the structure of the wavefunctions, 
we introduce the concept of projective symmetry group 
(PSG}~ which serves as the mathematical language to 
classify different symmetry protected FCI phases. 

With these theoretical preparations, in section Mil wc 
discuss one particular example, i.e. the checkerboard lat- 
tice mode l 11 ' 12 and write down two SU{m) FCI wave- 
functions and two Z m FCI wavefunctions in distinct uni- 
versality classes for v = 1/3. These wavefunctions sup- 
port the same <r xy — quantized Hall conductance 
and similar topological properties. They are character- 
ized by different PSGs in the bulk. These states can all 
serve as candidate states for the FCI state found in nu- 
merical simulation s 12 ' 13 ' 15 . Which state is realized in the 
simulated mode l 12 ' 13 ' 15 would be determined by energet- 
ics. Because our proposed wavefunctions has the form of 
a Slater determinant and can be effectively implemented 
by variational Monte Carlo approach, the energetics of 
the proposed states can be studied by future numerical 
investigation. In Appendix[G]we present another four ex- 
amples of distinct FCI phases in the honeycomb lattice 
model 9 -: two are SU(m) states and the other two are Z m 
states. We also propose spin-polarized FCI states with 
non-Abelian quasiparticlcs, which might be realized in 
nearly flat bands with Chern number C > 1. Such non- 
Abelian FCIs might be used to build a universal quantum 
compute r 29 ' 39 . 

In section lTvl we demonstrate that our parton construc- 
tion can be used to compute the topological ground state 
de generacy. This is particularly important for the Z m 
states, which belong to a new class of FQH wavefunc- 
tions. 

In section [V] we generalize our efforts to construct 
ground state wavefunctions of TRI FTIs. When the mix- 
ing between the up and down spins is weak in the elec- 
tronic hamiltonian, it is natural to generalize our spin 
polarized results to this case. For filling fraction v = ^- 
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(on average v = -j- for each spin), we present classes 
of SU{mY x SU (m)^ and x wavefunctions and 
discuss their properties including quasi-particlc statistics 
and ground state degeneracies. We also propose a new 
parton construction formalism which allows one to write 
down generic electron wavefunctions for TRI FTI states 
in the absence of spin conservation. We can deform such 
a generic TRI FTI wavefunction in the absence of spin 
conservation into a S^-conscrved TRI FTI wavefunction 
(where spin-t and spin- 1 decouple) by continuously tun- 
ing a parameter. Stability of such a state against pertur- 
bations are briefly discussed. 

II. SU(m) PARTON CONSTRUCTION OF 
SPIN-POLARIZED FRACTIONAL CHERN 
INSULATOR STATES 

A. A brief review of Laughlin's FQH state from 
SU(m) parton construction 

Soon after the experimental discovery of fractional 
quantum Hall (FQH) effects^, Laughlin proposed a se- 
ries of variational wavefunctions^ which were shown^ii 
numerically to be a very good description of FQH states 
at odd-denominator filling fraction v = 1/m. Later this 
idea of constructing trial wavefunctions was generalized 
to other filling fractions^—. An important lesson we can 
learn from Laughlin's wavefunction is as follows. With 
a fixed filling fraction (or a fixed number of flux quanta 
through the sample), the many-body wavefunction tends 
to vanish as fast as possible when two electrons approach 
each other so that the repulsive Coulomb energy between 
electrons could be minimized. As an example, Laugh- 
lin's state at v = 1/3 is nothing but the cube of the 
wavefunction for a filled lowest Landau level. We can 
construct this wavefunction by splitting an electron into 
three fermionic partons: 

c(r) = / 1 (r)/ 2 (r)/ 3 (r). (1) 

Naturally from (fTJ we can see each parton carries U(l) 
electric charge eo = e/3 where e stands for the electron 
charge. The electron wavefunction is obtained through 
the following projection 

N 

<M{rJ) - (0|n/i(r*)/2(r,)/ 3 (r 4 )|MF) (2 ) 

i=l 

where |0) represents the parton vacuum and \MF) can 
be any mean-field state of the three partons f 1,2,3- 
When each of the three partons occupy the lowest Lan- 
dau level (LLL) one immediately obtains the Laugh- 

lin's state $, =1/3 ({r 4 }) = [$Ui({rJ)] 3 = U t<J (^ " 

zj) 3 exp[— J2iLi \ z i\ 2 1 We have chosen the disc ge- 
ometry and the symmetric gauge. z\ = x% + \yi are com- 
plex coordinates, Ib = \Jhj\eB\ = l' B /y/3 is the electron 
magnetic length and l' B is the parton magnetic length. 



Since each kind of parton occupies a LLL, the electro- 
magnetic response of the FQH state $> v= i/ 3 ({ri}) = 

[$' ly=1 ({ri})] 3 is characterized by Hall conductivity 

Cxy = 3-(|) 2 //i = hjr- This reproduces the cor- 
rect filling fraction and the many-body Chern number. 
Note that electron operator c(r) in (fTJ is invariant un- 
der any local SU(3) transformation on the three par- 
tons (/i,/2,/3) T - The mean-field Hamiltonian density 
describing the Laughlin state is^ 

Umf = ^ E /«( r )( - iV - eoA(r)) 2 / Q (r) (3) 

where m* is the effective mass of each parton. This 
mean-field Hamiltonian preserves the SU(3) gauge sym- 
metry and partons will also couple to a SU(3) internal 
gauge field. Its effective theory is the SU(3)i Chern- 
Simons gauge theory, which explains the 3-fold topologi- 
cal ground state degeneracy on a torus^. These /t par- 
tons are nothing but charge e/3 quasiparticle excitations 5 
of Laughlin state. Indeed after projection ([2]) the three 
species of partons f 1.2,3 becomes indistinguishable thanks 
to the internal SU (3) symmetry: each / parton creates a 
charge —e/3 quasihole upon acting on the ground state 
\MF). It is straightforward to verify that the following 
wavefunction 

$e({ri}|w l!2)3 ) = (4) 
(0|/i(wi)/ a (w a )/ 3 (w 3 ) IltT 1 h{*i)M*i)h(?i)\MF) 

reproduces the Laughlin wavefunction with three quasi- 
holes at wi.2,3 up to a constant factor. Hence these 
partons are indeed charge e/3 anyons obeying fractional 
statistics with statistical angle 61/3 = -|. 

B. Z m FCI state and its quasiparticles from SU(m) 
parton construction 

Since the three seemingly-different partons /i,2, 3 are 
essentially the same quasihole excitations with the same 
quantum numbers, physically it is attempting to include 
the tunneling terms f\fp, a 7^ /? in the mean-field Hamil- 
tonian. By mixing different partons, these terms will 
break the internal SU (3) gauge symmetry down to a a 
subgroup of SU(3), which is Z 3 , the center of the SU(3) 
group, in the most generic case where fyfp, V a 7^ (3 
terms are present. In general the projected Z3 wavefunc- 
tion @ is different from its parent projected SU(3) wave- 
function. For a 2-D electron gas in a magnetic field, how- 
ever, people usually focus on the LLL within which the 
many-body wavefunction is an analytic function (e.g. in 
the symmetric gauge on a disc). It is straightforward 
to show that as long as the mixing terms act inside the 
Hilbert space of LLL, the corresponding electron wave- 
function for a Z3 state remains the same as that of 
its parent SU(3) state. This is because the parton wave- 
function describes a state with LLL fully filled. Mixing 
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between different partons within the LLL Hilbert space 
only gives a unitary transformation of basis and does not 
modify the parton wavefunction. For a lattice model, it 
is natural to consider mixing terms acting between all 
bands (rather than within the filled bands), and the cor- 
responding electron wavefunction of a Z% state will be 
a different wavefunction from that of its parent SU(S) 
state. For a filling fraction v = 1/m, our discussion 
straightforwardly generalizes to the corresponding Z m 
(the center of the SU(m) group) state and its parent 
SU(m) state. 

To our knowledge, the Z m parton states of FQHE have 
not been proposed before. For this new class of wave- 
functions, several natural questions need to be answered. 
What are the quasi-particles in the Z m state? What is 
the low-energy effective theory of the Z m state? Will it 
preserves the topological properties, such as ground state 
degeneracy? We answer these questions in this paper and 
find the topological properties of the Z m states are iden- 
tical to the SU(m) states. Their difference lies in the 
projective symmetry group, which is protected by lattice 
symmetry. In general, Z m states and SU (m) states both 
serve as candidate ground states for the FCI states of a 
v = 1/m filled band with Chern number one. 

To begin with, let us consider the quasiparticle exci- 
tations in a Z m state. The physical quasiparticle excita- 
tions in a Z m state are constructed by inserting fluxes in 
the mean- field ansatz of f\f a terms and simultaneously 
creating vortices (or defects) in the Higgs condensates 
(flfp), a ^ (3. In 2-D, because >K X {SU{m)/Z m ) = Z m , 
these defects arc the point-like vortices carrying Z m 
gauge fluxes. Because the Z m flux can be considered 
to be localized in a single plaqucttc, one can effectively 
interpret it as a overall U(l) gauge flux of all the /- 
partons. Namely when a /-parton winds around a fun- 
damental vortex, it experiences a 27r/m flux. Due to the 
Chern numbers of the filled parton bands, this vortex 
also binds with a single /-parton gauge charge and thus 
carries electric charge e/m. Because this object carries 
both flux and gauge charge, the fractional statistical an- 
gle 9 = n/m results. These are exactly the same electric 
charge and statistics that a quasiparticle carries in the 
SU(m) state. We conclude that the topological proper- 
ties of quasiparticles are the same in both the Z m and its 
SU(m) parent state. 

Following the above discussions, we can write down 
the wavefunctions with low-energy anyonic excitations in 
a Z m state. At filling fraction v = 1/m, in order to create 
one quasiparticle Q\ at wi and its antiparticle Qi at W2, 
we need to insert 2'nk/m flux in a plaqucttc P Wl at posi- 
tion wi and —2nk/m flux in a plaquette P W2 at position 
W2 with k = 1, • ■ ■ , to. Qi carries 2-nk/m flux and ek/m 
charge while Qi carries —2irk/m flux and — ek/m charge. 
Both Qi and Qi have statistical angle 8 = k 2 n/m and 
their mutual statistical angle is 0' = —k 2 n/m. They are 
realized by creating e l27rfc / m phase shift for all mean- field 
amplitudes on the string connecting two plaquettes P Wl 
and P W2 , on top of the mean-field ansatz for the ground 





FIG. 1: (color online) A pair of anyonic quasiparticle and 
its antiparticle in a Z m state on the honeycomb lattice. The 
dashed line denotes the string of e l27rk / m phase shift connect- 
ing the two plaquettes where quasiparticle Qi and its antipar- 
ticle Q2 are located. On top of the ground state mean-field 
ansatz, any mean-field bond crossing the string should pick 
up a phase shift e l27rfc / m Here we only demonstrate the phase 
shift of nearest-neighbor (NN) mean-field amplitudes by triple 
arrows in the figure. 

state. An example of such a pair of quasiparticle and 
its antiparticle in a Z m state on honeycomb lattice is 
schematically shown in FIG. [1] The corresponding elec- 
tron wavefunction is obtained by the projection of this 
new mean-field ansatz to the electronic degrees of free- 
dom. When to = 3, this projection is given by Eq.([2]). 

The ground state degeneracy of a Z m FCI state at 
v = 1/m on a torus can also be understood once we 
know its quasiparticle statistic a 31 ' 46 . Consider the fol- 
lowing tunneling process 71 : a pair of quasiparticle (with 
flux 2ix /m and charge e/m) and its antiparticle (with flux 
— 27r/m and charge —e/m) are created and the quasi- 
particle is dragged around the non-contractible loop X\ 
along x\ direction on the torus before it is finally an- 
nihilated with it anti-particle. This tunneling process 
will leave a string of e l27r / m phase shift (as shown in 
FIG.[T]) along this loop X\, therefore has the same phys- 
ical effects as adiabatically inserting a 2-k /to flux in the 
non-contractible loop X 2 along xi direction on the torus. 
Note that when the quasiparticle-anti-quasiparticle pair 
carries flux ±2Kk/m and charge ±ke/m the correspond- 
ing tunneling process is realized by 71 ■ Similarly we can 
define a tunneling process l~i by dragging the fundamen- 
tal quasiparticle around non-contractible loop Xi once, 
which is physically equivalent to inserting 2-k /m flux in 
non-contractible loop X±. In the thermodynamic limit, 
the Hilbert space of degenerate ground states should be 
expanded by these tunneling processes 4 ^. The two tun- 
neling operators satisfy the following "magnetic algebra" : 

71 T 2 = TiTie 2 ^" 1 (5) 
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This is straightforward to understand from the point of 
view of Aharonov-Bohm effect. Another way to under- 
stand it is because the tunneling process T1T2 
can de continuously deformed into two linked loops 3 -^ and 
corresponds to a phase of 29i/ m , where 0\/ m = ir/m is 
the statistical angle of the fundamental quasiparticle. All 
degenerate ground states can be labeled by e.g. eigenval- 
ues of unitary operators 71 and 7~™ (since they commute 
with each other) . In this basis T2 acts like a ladder oper- 
ator and changes the eigenvalue of 71 by a phase e l27r / m . 
In this way one can see the ground state degeneracy of 
a Z m state on torus is m-fold. We can easily generalize 
this discussion to a genius-*? Riemann surface with g pairs 
of non-contractible loops and the corresponding ground 
state degeneracy is m 9 -fold. This is consistent with the 
ground state degeneracy calculated from the low-energy 
effective theory as will be shown in section HVl in a formal 
way. 

Because the discussion on low-energy effective theory 
of the Z m state involves more technical details, we post- 
pone it to Section llVl where we compute its ground state 
degeneracy. We'll show that the ground state degener- 
acy of a Z m state is the same as that of a SU{m) state: 
m 9 -fold on a genus-g Riemann surface. 



C. Regarding lattice symmetries 

In the numerical simulations of FCI phase s 12 ' 13 , only 
the m-fold topological degeneracy of v = 1/m FCI is ob- 
served on torus. This indicates that the FCI wavefunc- 
tions respect the full lattice symmetry, since otherwise 
there should be extra degeneracies due to lattice symme- 
try breaking. This motivates us to write down the fully 
symmetric FCI wavefunctions. 

In the following we outline the general strategy to con- 
struct fully symmetric FCI states on a lattice in the par- 
ton approach. Here we focus on spin-polarized FCI states 
with filling fraction v = 1/m through the SU(m) parton 
construction. The electron operator is given by 



c(r)=n/«(r) 



(G) 



where r is the coordinate of a lattice site. For sim- 
plicity we assume there is only one orbital per lattice 
site. As mentioned earlier, this parton construction 
has a local SU (to) symmetry since the electron oper- 
ator c(r) is invariant under any local transformation 
U(r) -> E^G Q(3 (r)^(r) where G(r) e SU(m). A 
generic parton mean-field ansatz is written as 



^ = ^^/t(r)M afJ (r|r')/,(r') 



(7) 



a/3 



where M(r|r') = M^(r'\r) is a m x m matrix assum- 
ing there are one electron orbital per site. Under a lo- 
cal SU(m) gauge transformation {G(r)} it transforms as 



M(r|r') -» G(r)M(r|r')G t (r'). Again once we obtain a 
mean-field state \MF) with the right filling number from 
, the corresponding electron wavefunction is obtained 
through 



N 



N 



<M{rJ) = (0|nc(ri)|MF) = (0| \ \ / Q ( ri )|MF),(8) 



whose explicit form is a Slater determinant as given later 
in HQ]). 

Not all parton mean-field ansatzs correspond to v = 
1/m FCI states. Let's start from a SU(m) mean- 
field state with M aP (r\r') = S a ^T(r\r'): H™ F = 
E r ,r' fi{ v ) T a.A r \ T ')fa{*') where each flavor of the par- 
ton has the same filling number as the electron. For 
v = 1/m FCI states in topological flat bands, the fill- 
ing fraction is such that on average there is one elec- 
tron (hence one parton with each flavor) per m unit cells. 
If the mean-field ansatz has explicit lattice transla- 
tion symmetry, however, the corresponding state with 
v = 1/m filling would most likely be a gapless metallic 
stated since only a fraction (1/m) of the lowest band is 
filled. How to construct a gapped mean-field ansatz of 
FCI with filling fraction 1 / to? 

The answer lies in the SU(m) gauge structure of the 
parton construction © . This gauge structure allows the 
physical (lattice) symmetry to be realized projectively in 
the parton mean-field ansatz, which gives rise to a sym- 
metric electron wavefunction after projection. Briefly 
speaking, the mean-field state itself can explicitly break 
lattice symmetries (such as lattice translations) while the 
electron wavefunction after projection ((HJ) remains fully 
symmetric. 

By inserting e.g. 2-k /m flux in each original unit cell, 
one can enlarge the unit cell by m times. Therefore the 
corresponding mean-field state with filling v = 1/771 is 
a state filling the lowest m bands of mean-field ansatz 
([7]) . If each of the m lowest bands have a Chern number 
+1, the mean-field state filling these m bands would have 
total Chern number +m, and the corresponding Hall con- 
ductivity is 



>xy 



,1,2 ^ le* 
m-{ — ) ■ — = — — 

m h m h 



(9) 



because each parton carries U(l) electric charge e/m. 
This gives the correct electromagnetic response of a v — 
1/m spin-polarized FCI state. 

Here the mean-field Hamiltonian ([7]) explicitly breaks 
lattice translation symmetry due to the unit cell enlarge- 
ment, but as long as the translated mean-field ansatz can 
be transformed to the original ansatz by a local SU (m) 
gauge rotation, the corresponding electron wavefunction 
© still respect the translation symmetry. This is because 
any two mean- field ansatzs differ by a local SU(m) gauge 
rotation give exactly the same electron wavefunction af- 
ter projection. Similarly, even the mean-field Hamilto- 
nian breaks other lattice symmetries such as the point 
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group symmetry, the electron wavef unction after projec- 
tion still can be fully symmetric. The mean-field ansatz 
simply forms a projective representation of the symme- 
try group. The mathematical framework of construct- 
ing fully symmetric electron wavefunctions based on par- 
ton mean-field ansatzs is the projective symmetry group 
(PSG), which will be introduced shortly. The technique 
of enlarging the unit cell by m times without physically 
breaking any lattice symmetry will be generalized to the 
case of time-reversal-invariant FTI states with filling frac- 
tion v = 2/m in section fVl 

Following this strategy, we always require that the par- 
ton mean-field ansatz of v = 1 /m FCI state breaks lattice 
translation symmetry explicitly and enlarges the unit cell 
by m times, so that the resultant mean-field state is an 
insulator. The partons will fill the lowest m bands of 
the mean-field spectrum and the corresponding electron 
state after projection would still be gapped. Now that 
the number of momentum points of each band in the (re- 
duced) 1st Brillouin zone equals the electron number N, 
we can see the electron wavefunction ((5J) is nothing but 
a Slater determinant 



"M{r;})z m = detWmNxmN, 

W(a-l)N+i,(n-l)N+j = 0k"(r™)> 

a, n = 1, • • • , m; i,j = l,---,N. 



(10) 



where 0k™ (rf) represents the eigenvector of mean-field 
Hamiltonian ([7]). To be specific, 0k™ (rf) corresponds 
to the f a parton component in momcntum-kj single- 
particle eigenvector of the bottom-up n-th band. Here 
a, n = 1, • ■ ■ ,m and i,j = 1, ■ ■ • ,N where N is the to- 
tal electron number at filling fraction v = 1/m. Note 
that for a SU(m) mean-field ansatz (j2"5)) in the absence 
of mixing terms, the lowest m bands are all degenerate 
and we have 0k™ (r a ) = 0k(r)<S n , a . The corresponding 
electron wavefunction (I10p reduces to the product of m 
copies of a Slater determinant: 



$ e ({rJ)si/(m) = { 



dct 



0ki( r i)' 
0kj (r 2 )' 



fa) 
fa) 



.^(rjv)---^^). 



(ii) 



where 0k.,- (^i) is the momentum- kj single-particle eigen- 
vector of parton mean-field Hamiltonian H^ IF = 
E r , r ' / < £(r)T(r|r')/ a (r') with Va = 1, • • • , m. This is a 
lattice version of Laughlin's state in free spaced. However 
once we add lattice-symmctry-prcserving parton mixing 
terms which breaks gauge symmetry from SU(m) to Z m , 
the electron wavefunction (jTUJ) of a Z m FCI state, as well 
as its projective symmetry group which will be intro- 
duced shortly, will immediately become different from 
its parent SU(m) state (fTTj) . We emphasize again that 
only when the unit cell is enlarged by m times, we will 
have the same number of momentum points kj as the 
electron number N. The mean-field amplitudes can be 
determined by variational Monte Carlo study of the en- 
ergetics of electronic wavefunctions (fT0"|) . 



Considering flux insertion in order to enlarge the unit 
cell in the mean- field ansatz ([7)1. another question arises: 
can there be more than one way of inserting fluxes into 
plaquettes without breaking physical lattice symmetries? 
If yes, how to classify different mean-field ansatzs ([JJ? 
The answer of the 1st question is positive and to answer 
the 2nd question, we need to introduce a mathematical 
structure: projective symmetry group (PSG)^£ in order 
to characterize different "universality classes" of symmet- 
ric FCI states. PSG classifies different mean-field ansatz 
which forms a projective representation of the physical 
symmetry group. In the following we give a brief intro- 
duction of PSG. 

Note that there is a many-to-one correspondence be- 
tween parton mean-field states and physical electron 
states due to the above projection operation: any two 
parton mean-field states related to each other by a 
SU(m) gauge transformation {G(r)} correspond to the 
same electron state. As a result, although the physi- 
cal electron state preserves all lattice symmetry, its par- 
ton mean-field ansatz may or may not explicitly preserve 
these lattice symmetries. The physical lattice symmetries 
are realized projectivcly in the mean-field ansatz. More 
precisely, in a generic case the parton mean-field ansatz 
(UJ should be invariant under a combination of lattice 
symmetry operation U and a corresponding gauge trans- 
formation {G[/(r) e SU(m)}: 

M(U(r)\U(r')) = G c/ (C/(r))A/(r|r')G^(C/(r')) (12) 

Different universality classes of parton mean-field ansatzs 
arc characterized by different PSGs2£, i.e. different 
SU(m) gauge transformations {Gu} associated with 
symmetry operations U: 

PSG = {Gu(r)U\U e symmetry group} (13) 

The low-energy gauge fluctuation of a mean-field ansatz 
is controlled by its invariant gauge group2£ (IGG) 

IGG = {G e g SU{m)\G e M(r\r')Gl = M(r|r'),V r,r'} 

where e represents the identity operator of the (lattice) 
symmetry group (SG). In other words, IGG is a subgroup 
of the internal gauge group (which is SU(m) here) that 
keeps the mean-field ansatz ([7J invariant. Hereafter we 
would call a parton mean-field state with e.g. IGG = 
SU(m) state a SU(m) state. We can see that the IGG of 
mean-field ansatz (UJ always contains the following Z m 
group as a subgroup: 



Z n 



■In 



\a = l,2, 



(14) 



where I m xm is the mxm identity matrix. This Z m group 
is the center of the SU(m) group. A mean-field ansatz 
{M(r|r')} with IGG = Z m is called a Z m stated. The 
low-energy theory of this Z m state will be described by 
fcrmionic partons f a interacting with Z m gauge fields. 

The classification of PSGs with IGG = SU{m) (which 
we call SU (m) PSGs in this paper) are easy to carry out. 



FIG. 2: (color online) The lattice structure and symmetry 
group generators of (a) honeycomb lattice and (b) checker- 
board lattice. Si and S2 are primitive lattice vectors. Each 
site is labeled by coordinate (x,y,s): f = xa\ + ya,2 corre- 
sponds to the position vector of its unit cell and s — 0/1 are 
the sublattice indices for A/B sublattices. Ce represents n/3 
rotations along z-axis around the honeycomb plaquette cen- 
ter, er and P x , P y , P xy are n rotations along the axis plotted 
in the figure. 



The only gauge invariant quantities of a SU (m) ansatz 
is the gauge-invariant flux through each plaquette, which 
must belong to the center of the SU(m) gauge group, 
namely the Z m group in Eqll4[ because otherwise the 
SU(m) gauge group would be broken and IGG cannot 
be SU(m). Two SU(m) states have the same PSG if and 
only if they have the same Z m gauge flux in each given 
plaquette. Therefore distinct SU(m) PSGs have different 
Z m gauge flux pattern and vice versa. 

The classification of PSGs with IGG = Z m (which we 
call Z m PSGs in this paper) involves more technical de- 
tails and we leave this analysis for the honeycomb lattice 
models and the checkerboard lattice mode l 11 ' 12 in Ap- 
pendix [B] and [D] 



III. EXAMPLES 

The first part of this section shows four concrete exam- 
ples of mean-field ansatzs corresponding to spin-polarized 
FCI states with filling fraction v = 1/3 in the checker- 
board lattice mode l 11 ' 12 . These include two Z3 ansatzs 
and their two SU(3) parent states. As discussed in the 
previous section and proved in Appendix lATlDl although 
the mean-field ansatz explicitly breaks lattice translation 
symmetry by tripling the unit cell, the physical elec- 
tron state after projection ([8]) preserves all lattice sym- 
metry. Four examples for the honeycomb lattice model 
with filling fraction v = 1/3 are displayed in Appendix 
[Gl We show the Hall conductance of all these states are 
&xy = e 2 /(3h). 

In the second part of this section, we present a sce- 
nario which might realize non-Abelian FCI with IGG = 
SU{m) by partially filling nearly flat bands with Chern 
number C > 1, and two examples of such states in the 
SU(3) parton construction are presented in Appendix [HI 



FIG. 3: (color online) The 1st Brillouin zone of (a) honey- 
comb lattice and (b) checkerboard lattice, bi and 62 denoted 
by arrows are reciprocal lattice vectors satisfying bi ■ Sj — Sij 
where 01,2 are primitive vectors shown in FIG. [2] The solid 
line (red) denotes the original 1st Brillouin zone (BZ) of the 
two lattices while bold dashed line (blue) denotes the re- 
duced 1st BZ when ±27r/3 fluxes are inserted in each unit 
cell. The area of the reduced BZ is only one third of the orig- 
inal BZ. (a) For the honeycomb lattice, the three circles have 
momenta F : (fa, fa) = (0,0), K : (fa, fa) = (77r/6,7r/3) 
and M : (fa, fa) = (tt,0), where k = fa&i + fab2- (b) In 
the case of checkerboard lattice, the three circles have mo- 
menta T : (fa, fa) = (0,0), K : (fa, fa) = (7r/3, 7r) and 
M : (fa, fa) = (0, 7r), where k = fa&i + fa&2. 



A. Four examples of spin-polarized FCI states with 
v = 1/3 in the checkerboard lattice model 

The checkerboard lattice mode l 11 ' 12 has been shown 
to support nearly flat bands with non-zero Chern num- 
bers. Its symmetry group is shown in Appendix [C] and 
FIG.[2j Each lattice site is labeled by coordinate (x, y, s) 
as shown in Appendix [Cj Recently there are numerical 
evidenc o 12 ' 13 ' 15 of FCI states with filling fraction v = 1/3 
and v = 1/5 in this model. By SU(m) parton con- 
struction, we use PSG to classify different Z m mean-field 
ansatz (m is an odd integer) as shown in Appendix [Dl In 
the following we show two Z m mean-field states that be- 
long to different universality classes. They correspond to 
two gauge-inequivalent solutions of (|D13[) when m = 3. 
We show mean-field amplitudes up to the next next near- 
est neighbor. Their two parent SU(3) states also have dif- 
ferent PSGs because they have different patterns of the 
SU(3) gauge invariant fluxes. All four states are candi- 
dates of the FCI state found in the exact diagonalizations 
of checkerboard lattice mode l 12 ' 13 ' 15 . 



1. The Z3 FCI state CB1 with a xy = 1/3 and its parent 
SU(3) state 

In the Z3 FCI state CB1 the gauge transformations 
Gu(x,y,s) associated with lattice symmetries U are 
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FIG. 4: (color online) Mean-field ansatz of the parent SU(3) 
parton states associated with FCI state CB1 on checkerboard 
lattice. The solid lines on NN bonds represent complex hop- 
ping amplitude a along the direction of the arrow. Solid lines 
on NNN bonds represents real hopping amplitude f3 x . The 
dashed lines on NNN bonds represent complex hopping am- 
plitude e'^^fly along the direction of the arrow. The triple 
arrow means the original hopping amplitude along its direc- 
tion should be multiplied by a phase factor 7712. We only show 
up to NNN terms and all NNNN terms can be obtained from 
Note that in the mean-field ansatz the lattice transla- 
tion along Si direction is explicitly broken by flux insertion 
and the unit cell is tripled. However the lattice translation 
symmetry is preserved in the corresponding electron states 
after projection iJHJ. 



CB1: parent SU(3) state 
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FIG. 5: (color online) The parton band structure of the par- 
ent 5(7(3) ansatz ^ of CB1 state. Each band is 3-fold 
degenerate, corresponding to 3 parton flavors /i,2,3. We plot 
the dispersion along F —¥ K — > M — > F as shown in FIG. 
|3jb). Hopping parameters are chosen as a = e"^ 12 for 
NN, p x = -0.2 = p y for NNN and 7 = 0.1 for NNNN. The 
Chern numbers of the 6 bands are {1, — 2, —2, 4, —2, 1} in a 
bottom-up order. 



listed below: 

Gt 2 (x, y, s) = 7 3x3 , G Tl (x,y,s) = r)\ 2 I Sx3 , (15) 

Gp m (x,y,s) = r){%I 3 x3, 

Gp y (x, y, s) = 73x3, G PxB (x, y, s) = ??if I 3 x3- 

where 7712 = exp(— i27r/3). T\,i are lattice translation 
operations and P x _ y _ xy arc reflections. They are schemat- 
ically shown in FigJ5Jb) and defined mathematically in 
Appendix [D] As shown in Appendix [D] the symmetry al- 
lowed mean-field amplitudes are: 

(I) For nearest neighbor (NN) amplitude u a = 
M(0,0,1|0,0,0) 

u a = (16) 

i.e. u a can be any complex symmetric 3x3 matrix. All 
other NN amplitudes can be generated from u a by sym- 
metry operations through (|12p . 

There are two independent NNN amplitudes: 
(Ha) For next nearest neighbor (NNN) amplitude 
u^ = M(l,0,0|0,0,0) 

u Px = u T px = u* 0x (17) 

i.e. up x can be any real symmetric 3x3 matrix. Half of 
NNN mean- field amplitudes can be generated from up x 
by symmetry operations through ()12p . 

(lib) For next nearest neighbor (NNN) amplitude 
^ = 1/(0,1,010,0,0) 

u P y = e i * /3 u y (18) 

where up v can be any real symmetric 3x3 matrix. Half 
of NNN mean-field amplitudes can be generated from up y 
by symmetry operations through (|12p . 

(Ill) For next next nearest neighbor (NNNN) ampli- 
tude up = A/(1,1,0|0,0,0) 

u 7 = u\ (19) 

i.e. u 7 can be any Hermitian 3x3 matrix. All other 
NNNN mean-field amplitudes can be generated from w 7 
by symmetry operations through \12\ . 

These mean- field amplitudes of u a , up x , up y and w 7 
should be treated as variational parameters. Their opti- 
mal values which minimize the variational energy of elec- 
tron wavefunction (|10p can be determined in variational 
Monte Carlo simulations. 

The corresponding parent SU (3) state has u a = 
ahx3, u ax = p x I 3 x3, u ay = e 17T/3 /3 y I 3x3 and u 7 = 
7^3x3- Choosing parameters a = e 17r / 12 , (3 X = —0.2 = @ y 
and 7 = 0.1 we have the Chern numbers of the 6 bands 
for each parton species as {1, — 2, — 2, 4, — 2, 1} and the 
lowest band is well separated from other bands, as shown 
in FIG. 03 This qualitative band structure persists for a 
large parameter range. Each band of the SU (3) parton 
ansatz is 3-fold degenerate, corresponding to the 3 parton 
flavors /i,2,3- By adding small parton mixing terms to 
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CBI: Z3 state 



5 







-5 



; 


— - — 






















^ i " =T 


3 










r k 


M 


r 



FIG. 6: (color online) The parton band structure of a Z3 
mean-field ansatz (J7J): the CBI state. We plot the dispersion 
along r -> K -> 71/ T as shown in FIG. E^b). Note that 
each 3-fold degenerate band in its parent 5(7(3) ansatz (see 
FIG. [5} now splits into 3 non-degenerate parton bands. All 
the 3 lowest bands have Chern number +1. 

the SU (3) mean-field state, each 3-fold degenerate band 
splits into 3 non-degenerate parton bands in a Z 3 mean- 
field state, as shown in FIG. [51 By filling the resulting 
3 lowest bands (all with Chern number +1) we obtain a 
Z 3 FCI state whose Hall conductivity is a xy = 1/3 in the 
unit of e 2 /h. 

2. The Zz FCI state CB2 with a xy = 1/3 and its parent 
5(7(3) state 

In the Z 3 FCI state CB2 the gauge transformations 
Gu{x, y, s) associated with lattice symmetry U are listed 
below: 

Gt 2 (x, y, s) = 7 3x3 , G Tl {x,y,s) = Tf 12 I 3x3 , (20) 

Gp x (z, y, s) = ?^2~ 1);r/ 3x3, 
Gp y (x,y,s) = ri\ 2 I 3x3 , G P:ry (x,y,s) = ?7if/ 3 x3- 

where 7712 = exp(— i27r/3). As shown in Appendix [D] the 
symmetry allowed mean-field amplitudes are: 

(I) For nearest neighbor (NN) amplitude u a = 
M(0,0, 1|0, 0, 0) 

u a = u T a (21) 

i.e. u a can be any complex symmetric 3x3 matrix. All 
other NN amplitudes can be generated from u a by sym- 
metry operations through (jT2J) . 

There are two independent NNN amplitudes: 
(Ha) For next nearest neighbor (NNN) amplitude 
u /jB = Af(l ) 0,0|0,0 J 0) 

u fjx = e^^u 0x (22) 

where up x can be any real symmetric 3x3 matrix. Half 
of NNN mean-field amplitudes can be generated from up x 
by symmetry operations through (|12[) . 
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FIG. 7: (color online) Mean-field ansatz of the parent S?7(3) 
parton states associated with FCI states CB2 and CB3 on 
checkerboard lattice. The solid lines on NN bonds represent 
complex hopping amplitude a along the direction of the ar- 
row. Dashed lines on NNN bonds represents real hopping 
amplitude f3 y . Solid lines on NNN bonds represent complex 
hopping amplitude e ln ' 3 /3 x along the direction of the arrow. 
The triple arrow means the original hopping amplitude along 
its direction should be multiplied by a phase factor 7712. We 
only show up to NNN terms and all NNNN terms can be 
obtained from (JT2J ■ Note that in the mean-field ansatz the 
lattice translation along Si direction is explicitly broken by 
flux insertion and the unit cell is tripled. However the lat- 
tice translation symmetry is preserved in the corresponding 
electron states after projection ([8|. 

(lib) For next nearest neighbor (NNN) amplitude 
U0„ = M(O,1,O|O,O,O) 

up y = u T py = u} y (23) 

i.e. upy can be any real symmetric 3x3 matrix. Half of 
NNN mean-field amplitudes can be generated from up y 
by symmetry operations through (|12p . 

(Ill) For next next nearest neighbor (NNNN) ampli- 
tude up = A/(1,1,0|0,0,0) 

u-y = u\ (24) 

i.e. « 7 can be any Hermitian 3x3 matrix. All other 
NNNN mean-field amplitudes can be generated from m 7 
by symmetry operations through p2[) . 

The corresponding parent SU (3) state has u a = 
ahx3, u ax = e 17r / 3 /3 x I 3x3l u ay = /3 y I 3x3 and u 7 = 
7/3x3. From FIG. [7] we can see that the pattern of fluxes 
in CB2 state is different from that in CBI state. For 
example, in CBI state — 27r/3 flux (shown by triple ar- 
rows) are inserted in one half of the NN square plaquettes 
(those enclosed by two horizontal blue sites and two ver- 
tical red sites) while in CB2 state — 27r/3 flux are inserted 
in the other half of the NN square plaquettes (those en- 
closed by two horizontal red sites and two vertical blue 
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CB2: parent SU(3) state 
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FIG. 8: (color online) The parton band structure of the par- 
ent S(7(3) ansatz of CB2 state. Each band is 3-fold 
degenerate, corresponding to 3 parton flavors /i,2,3- We plot 
the dispersion along Y — > K — > M — > V as shown in FIG. 
|3fb). Hopping parameters are chosen as a = e lv ' 12 for 
NN, p x = 0.1 = P v for NNN and 7 = 0.09 for NNNN. 
The Chern numbers of the 6 bands are {1, 1, 1, —5, 1, 1} in 
a bottom-up order. 

sites). Notice that these two types of NN square plaque- 
ttes are inequivalent, not related to each other by any 
lattice symmetry (such as translation or mirror reflec- 
tion) in the checkerboard lattice modelii. Therefore CB1 
and CB2 are different mean-field Z3 FCI states belong- 
ing to distinct universality classes. Choosing parameters 
a = e i7r / 12 , f} x = 0.1 = P y and 7 = 0.09 we have the 
Chern numbers of the 6 bands for each parton species as 
{1,1,1,— 5, 1,1}, and the lowest band is well separated 
from other bands, as illustrated in FIG. |3 This qualita- 
tive band structure persists for a large parameter range. 
Each band of a SU(3) parton ansatz is 3-fold degenerate, 
corresponding to the 3 parton flavors j 1,2,3 ■ By adding 
small parton mixing terms to the SU(3) mean-field state, 
each 3-fold degenerate band splits into 3 non-degenerate 
parton bands in a Z3 mean-field state. By filling the re- 
sultant 3 lowest bands (all with Chern number +1) we 
obtain a Z3 FCI state whose Hall conductivity is again 
°xy = 1/3 in unit of e 2 /h. 



B. Possible non-Abelian states by partially filling a 
nearly flat band with Chern number C > 1 

It has been show n 45 ' 47 that in SU (m) parton construc- 
tion, when each parton species fills n Landau levels, the 
effective theory of the corresponding electron state is the 
SU(m) n Chcrn-Simons theory and the system has non- 
Abelian quasiparticle excitations when n > 1. Moreover, 
the non-Abelian quasiparticles of this state can be used 
as topologically-protected qubits in a universal quantum 
computer— as long as m > 2. This motivates us to pro- 
pose possible realization of non-Abelian FCI states re- 
alized in a partially-filled nearly flat band with Chern 
number C > 1. When this band is partially filled with a 



filling fraction e.g. v — 1/m, the Hall conductance of the 

corresponding FCI state would be a = ~ %- . If the lowest 
to parton bands all have Chern number C instead of +1, 
the SU(m) FCI state obtained by filling to lowest parton 

bands indeed has Hall conductance a — m^p 1 = — 4r 

h m h 

with en = e/m. This SU(m) FCI state can be a promis- 
ing candidate for v = 1/m-filled nearly flat band with 
C > 1. 

We discuss the C = 2 case as an example. In Ap- 
pendix [H] using the SU(3) parton construction, we show 
two examples of SU(3) FCIs, one on the honeycomb lat- 
tice and the other on the checkerboard lattice, both have 
3 lowest parton bands with Chern number +2. These 
two SU(3) FCIs with v = 1/3 and a = are non- 
Abelian FCIs. Their low-energy effective theory of these 
SU(3) states is SU(3)2 Chern-Simons theory^ 5 -, featured 
by 6-fold ground state degeneracy on the torus and non- 
Abelian quasiparticle excitations^. These results indi- 
cate that once a nearly flat band with Chern number 
C > 1 is found, by partially filling it one may realize 
non-Abelian FCIs, which have the potential to build a 
universal quantum compute r 29 ' 39 . 



IV. EFFECTIVE THEORY AND GROUND 
STATE DEGENERACY OF SPIN-POLARIZED Z m 
FCI STATES 

As mentioned earlier, the partons in our construction 
not only couples to the external electromagnetic gauge 
field, but also couples to a dynamical internal gauge field, 
which is a SU(m) (Z m ) gauge field for a SU(m)(Z m ) 
state. The low energy effective theories of the partons 
coupled with the internal gauge fields control all the topo- 
logical properties of the systems. The topological prop- 
erties of a SU(m) state has been studied befor e) 44 ' 45 . In 
this section we will analyze the low-energy effective the- 
ory of spin-polarized Z m FCI state (|7|) from the SU(m) 
parton construction. To our knowledge, the Z rn FQH 
states have not been proposed and studied before. We'll 
try to answer the following questions: what is the ground 
state degeneracy of the Z m FCI state? Is it the same as 
or different from that of a SU(m) FCI state? 

We start from a SU(m) mean-field state which has 
been shown^ 5 - to describe the v = 1 /to Laughlin state in 
the continuum limit. Its mean-field ansatz is 

M a , p (r\T f ) = 6 a , p r(r|r') (25) 

In other words, there is no hopping between partons of 
different species and the to species of partons have ex- 
actly the same band structure. As shown in (jllj) its 
electron wavefunction is a lattice version of Laughlin 
state 5 -. Apparently this mean-field state doesn't break 
the SU (to) gauge symmetry which leaves the electron 
operator ((6|) invariant. Since here the partons couple to 
both the U(l) electromagnetic field and the SU (to) gauge 
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field, the Lagrangian writes 



J SU(rn) 



(26) 



/ dt{ Ea E r i/i (r, Wa(r, t) - Ea.r.H ^Ir') 

. e -ieo / r r ' A(x,t)-dxp [ e -i / r r ' 3(x,t)-dij /^(r'^) j ) 



where V means path-ordered integral. and a M are the 
[7(1) and the SU(m) gauge fields respectively. Strictly 
speaking they are both defined on the link of a lattice 
here. To linear order the above action can be written as 



L 



SU{m) 



(27) 



/ d*{ E Q E r /l(r, t) [id t f a (r, t) - Er' T(r\v')f a (v', t) 
-eoE r ^ (1) M)^(r,i) 
-Er(4 U{m) ) a ,,(r,t) <>,*)} 

where A^(r) stands for electromagnetic U(l) gauge field 
while a M (r) represents the internal SU(m) gauge field, 
eo = e/m is the electric charge of each parton. Here 
J u W(r) and J^ U(m] (r) are conserved 17(1) and SU(m) 
parton currents respectively. To be precise, Jo = 

Ea /a/a an d ( ) a p = flfp ~ S a> p E 7 tyf-y/™" In 

the long- wavelength limit the spatial components of the 
parton currents in momentum space (with momentum q) 
writes: 



(7(1) 
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Ek^k7kEa /l,k-q/2/ Q , k +q/ 2 ' 



(J^q (? ™ ) a ,p - Ek^k^k/lk-q^/s, 



q / 2 //3,k+q/2 



$a., 



jU{1) 



Since the partons form a band insulator, the band gap 
allows us to safely integrate out the partons and obtain an 
effective action L [A, a] for the gauge fields. Let's assume 
all the filled m lowest parton bands have Chcrn number 
+1. Upon integrating out partons {f a , /t}, the effective 
Lagrangian density writes 



(28) 



the first term corresponds to the quantized Hall con- 
ductance <T xy — me^/h. while the second term, i.e. a 
SU(m)i Chcrn-Simons term describes the low-energy 
gauge fluctuations. As shown in Appendix [El the Chern- 
Simons theory of SU(m) gauge field a M can be reduced 
to Chern-Simons theory of U(l) gauge fields af, I = 
1, • • • , m — 1. The gauge field configuration is given by 
( a n)a.p = Ejwl a \i9a p wnere 9 1 are m x m matrices 
defined in (|E2[) . In the ao — gauge, a\ and ai are con- 
jugate variables since the Lagrangian density for internal 
gauge fields a M writes 



_^ m — 1 

£sU(rn) 1 [°] = J" X] 7 ( J + ^ ( a l 9 '°2 _ a 2^a( 
/=1 



According to uncertainty principle, a{ and a\ cannot be 
determined simultaneously and we choose to fix the con- 
figuration of a\. Aside from these Z7 (1) gauge symme- 
tries, there are also discrete symmetries associated with 
essentially all permutations between partons (for details 
see Appendix |E| . Taking all these into account we can 
obtain the ground state degeneracy as the number of 
gauge- inequivalent configurations of {a^}. As shown in 
Appendix [El the m-fold degenerate ground states corre- 
spond to the following gauge field configurations: 



a 



0, (J =1,2, 

m— 1 2JLJL I 

2 — L 2 m> " 



1. 



-2); 

, m. 



(29) 



Physically this means once we insert 2irk/m flux in the 
hole along x\ direction of the torus for each parton, the 
original ground state is transformed into a different de- 
generate ground state. This is a "small" gauge transfor- 
mation for the partons since they transform as 



fa 



exp 



m—i 

i=l,2 



4dat 



(30) 



1=1 



Now we add Higgs terms M a ^(r\r') which break the 
original SU(m) gauge symmetry down to Z m . Does 
the corresponding Z m state have the same ground state 
degeneracy as a SU(m) state? The answer is posi- 
tive. In the long-wavelength limit we introduce the 
Higgs fields Q( g which carry no electric U(l) charge but 
carry the internal gauge charge. As an example, the 
/|(r)Mi jm (r|r')/ m (r') terms in the lattice model will in- 
troduce Higgs field 4>i. m (xi,X2) in the long-wavelength 
limit. The Higgs field </>i, m carries a 1 charge +1 for 
I = 1, • • • , m — 2 and a m ~ 1 charge +m. Likewise, for 
example, (j> 2 , m carries a 1 charge —1, a 1 charge +1 for 
I = 2, • • ■ , 771 — 2 and a m_1 charge +m. In general for 



a Higgs field 



J P,a 



a < (3 associated with mixing 



term f\fp has a 1 ^ charge Q T a g where 



Qi 







0, 

a — 1, 
+1, 

P, 

In the end one can see that 
can be viewed as adding a 
gauge field configurations 
action (|2"9")l . To be specific, 
with the presence of Higgs 
density for internal gauge 



or I > 



(31) 



I < a-2 
I = a-l 

a<I </3-2 
1 = 13-1 

the condensation of Higgs field 
potential in the phase space of 
to the SU(m)i Chern-Simons 
once we integrate out partons 
fields the effective Lagrangian 
fields becomes 

p] = £,SE/(m)i [ a ] 

iEi Qi^D^a,? 



Note that the above action is invariant under the follow- 
ing "large" gauge transformations 

(4,4)^(4 + ^,4 + ^1 
<t> a ,e [i27r Ei QiA 2 ^ + v -fr) 
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where pi^ are integers so that <p a ,i3 is a single-valued 
function on the torus. Besides there are other large gauge 
transformations as listed in (|E8[) . And all the discrete 
gauge transformations associated with permutations be- 
tween partons are also present, such as Pi 2 



h 



Upon integrating out the fluctuations of Higgs 
fields Sipa^p around their mean-field values <$> a ,fs m 
Ceffla 1 ,(p a .p} we have 



£z m [a] = C SU (m)i H - V[a[, a[] 



(32) 



The exact shape of potential U[a{,a2] depends on 
e.g. magnitudes of Higgs fields </> a ,/3; but it has certain 
robust features determined by the gauge charges ^ 
of Higgs fields 4> a ,p which condense^. More precisely, 
this potential is periodic in configuration space, with 
periodicity 2ir/Li for a{, I < m — 2 and 2ir/Li2 for 
2. It also has periodicity 27r/(mLi) for 



a\ , I < rn 



and l-KjiraL^) for a r 2 



The minima of this po- 



tential sits exactly on the configurations shown in (|2l?l) of 
the m-fold degenerate ground states of SU(m)i Chern- 
Simons theory. Besides these features associated with 
large gauge transformations, the potential V[a{, a?] are 
not invariant under the discrete symmetries associated 
with parton permutations such as Pi. 2- This is essen- 
tially because the introduced mixing terms (or Higgs con- 
densation) breaks the SU{m) gauge symmetry. The La- 
grangian (f3"2"j) actually describes the motion of particles 
in a magnetic fiel d 45 ' 49 and a periodic potential: the 7-th 
particle associated with a 1 ^ experiences a magnetic field 
of 7(7+1) flux quanta piercing through the torus. Due to 
the periodicity of potential Vfaf,^], the m-fold ground 
state degeneracy (as calculated in Appendix [Ej is still 
present when the gauge symmetry is broken from SU (m) 
down to Z m by introducing mixing terms between differ- 
ent partons (or Higgs fields). 

This can be understood physically: by threading a 
2-nk/m flux of gauge field a™ -1 in the hole along x\ di- 
rection on the torus, one creates a vortex (or 2ix phase 
winding) in the <t>i, m , I = 1, • ■ ■ ,m — 1 condensates in 
the non-contractiblc loop along the xi direction. This op- 
eration exactly corresponds to the tunneling process T 2 k 
mentioned in section [H] and will cost zero energy in the 
thermodynamic limit. Therefore the presence of Higgs 
fields will not lift the m-fold ground state degeneracy in 
the thermodynamic limit. 

We've shown that the ground state degeneracy of a 
Z m FCI state is m on a torus. The above analysis for 
the torus case can be easily generalized to study the 
ground state degeneracy on a genus-g Riemann surface. 
There are g pairs of non-contractible loops {A a ,B a \a = 
1, • • • , g} on a genus-g Riemann surface where each pair is 
just like the two non-contractible loops on a torus. Thus 
one can straightforwardly show that the ground state de- 
generacy of a Z m FCI state is m 9 on a genus-g Riemann 
surface. 



V. PARTON CONSTRUCTION OF 
TIME-REVERSAL-INVARIANT FTI STATES 

A. SU(my x SU(m) i ' parton construction of TRI 
FTI states 

We have focused on spin-polarized FCI states. Taking 
into account spin degrees of freedom, nearly flat bands 
with time-reversal-invariant (TRI) Z2 index-ii can exist. 
When a pair of bands carrying Z2 index are partially 
filled, can a fractionalized topological phase preserving 
both time reversal and lattice symmetries be realized? 
In principle the answer is yes. As a direct generaliza- 
tion of spin-polarized SU (m) and Z m FCI states, when 
the pair of nearly flat Z2 bands are filled partially with 
v = 2/m, by SU{m)^ x SU(m)^ parton approach we 
can construct a TRI fractionalized phase which we term 
as SU{ m y x SU(m)± and Z} n x Z^ FTIs. In a sim- 
ple way, a SU(m)^ x SU(m)± {Zj n x Z^) FTI wave- 
function with v = 2/m is a direct product of a spin- 
polarized SU{m) (Z m ) FCI state with a xy = e 2 /(mh) for 
spin f and its time-reversal counterpart: a spin-polarized 
SU{m) (Z m ) FCI state with a xy = —e 2 /(mh) for spin 
4,. Consequently in these FTI wavef unctions, entangle- 
ment between electrons with opposite spins is absent. We 
term these wavefunctions as spin-conserved FTI wave- 
functions, and will construct FTI wavefunctions involv- 
ing entanglement with opposite spins in the next subsec- 
tion. But as a fully gapped topological phase, this phase 
is stable at least when the mixing between spin j" and 
spin 4, is weak in the electronic Hamiltonian. As in the 
case of spin-polarized SU (m) and Z m FCI states, we still 
use the technique of enlarging the unit cell by m times, 
to guarantee that the ground state with the correct filling 
fraction is an insulator with a band gap. 

The mean-field ansatz of a generic spin-conserved FTI 
wavefunction is written as 

- Er^E^' fLAr)M a Ar,<y\r',a')f f3Ml (r'), 



M a ^(r|r')<5 CT)t + M*^(r|r')^ 



(33) 



where a, a' =t,4- are the spin indices. M a ^(r\r') can 
be any mean-field ansatz of a SU (m) (or Z m ) FCI state 
constructed in sectionlHland demonstrated in section Hill 
And the corresponding FTI state is a SU (m)^ x SU (m)^ 
(Z^ x Z^ n ) FTI. Apparently (|33)) is invariant under time- 
reversal transformation while preserving all the lattice 
symmetries. 

Again the iV-electron wavefunction (with N/2 elec- 
trons for each spin here in the 5' z -conserved case) is 
obtained by projection on mean-field state \MF) = 
\MFf) g> IMF,.) 



<I> e ({rJ};{rH) 



N/2 



N/2 

rW 1 



■)\MF), (34) 
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which is simply a product of two Slater determinants (|10p 
for spin up and down partons. Following our discussion in 
section Hvl the ground state degeneracy of this SU(m)^ x 
SU(m)^ (Z^ x Z^) FTI is m 2 on a torus and there are 
quasiparticle excitations (with both spin) with electric 
charge ±e/m. 

The SU(m)T x SU(m)± (or Zj n x Z^ n ) FTI wavefunc- 
tions constructed above explicitly conserve the both the 
t and J. electrons. However, in a spin-orbit coupled sys- 
tem where the S z conservation is not a symmetry, the 
true ground state wavefunction must involve mixings be- 
tween the two spin species. Of course it is possible that 
this true ground state is in the same universality class as 
those spin-conserved FTI states, because they are gapped 
stable topological phases. Nevertheless it is still interest- 
ing to explicitly write down a FTI state without spin- 
conservation. 

There is a natural question that needs to be an- 
swered in the current formalism: when spin mixing terms 
fa are present in the parton mean-field Hamilto- 

nian Eq |33l while preserving the time-reversal symmetry, 
is the corresponding electronic state EqJM] a TBI FTI 
wavefunction without spin conservation? The answer is 
negative. Below we study the properties of this state in 
details. 

The mean-field Hamiltonian including spin mixings be- 
tween partons is 

h mf = h mf + ^ [/t )t ( r )M a ,0(r,t |rU)/^u(r') 
+ h.c] (35) 

Note that upon mixing partons with different spins, 
e.g. for a Z^ x Z^ FTI the internal gauge symmetry is 
further broken from Z^ x Z^ (one Z^ for spin-f parton 
and another independent Z^ for spin-J, parton) 

/o,<r(r) -> e 1 ^ 2 'fa,a(r), Of, 04. = 1,2, m. (36) 

to a single Z m symmetry (for partons with both spin j" 
and 4)64: 

/ Q , CT (r) ^ e'^/^r), o = 1,2,- ■ • ,m. (37) 

In this paper, we define the FTI phase as the phase 
of matter which is in the same universality class as the 
direct product of two FCI states of opposite spins while 
preserving the time-reversal symmetry. Following this 
definition, the Z m state defined in Eg 1351 with spin mix- 
ings is not a FTI state. 

One clear difference between the Z m state Eq[3S] and 
the x Z^ FTI state is their quasiparticle statistics. 
In a x Z^ FTI state, the Z^ L fluxes and the Z^ fluxes 
bind with parton charges due to the Chern numbers of 
the parton bands, and are anyons with statistical angles 
= ±jjj. But in a Z m state Eol35l the Z m fluxes do not 
bind with parton charges, and are bosons. This indicates 
that the Z m state Eq 1331 which is described by a regular 
Z m lattice gauge dynamics with bosonic flux excitations, 



must separate from the Z^ m x Zj^ FTI state by a phase 
transition. 

Another observation that we have is that, in the 
SU(m)^ x SU(m)^ parton construction, the Z m state 
Eg 1351 cannot preserves both time- reversal symmetry and 
lattice translation symmetry simultaneously. Essentially 
the technique of enlarging unit cell m times by inserting 
fluxes fails to generate a translation symmetric wavefunc- 
tion when partons with opposite spins are mixed. In a 
TRI Z^ x Z^ FTI state, when 27r/m fluxes are inserted 
in each unit cell for spin 4" partons, an opposite —2ir/m 
fluxes arc inserted in each unit cell for spin 4- parton to 
preserve the time-reversal symmetry. After the partons 
with opposite spins are mixed, the gauge flux pattern of 
the partons no longer enjoys well-defined ±27r/m value 
per plaquette. As a result physical translation symme- 
try is broken. This simple argument dictates that, us- 
ing SU(m)^ x SU(m)^ parton construction, when par- 
tons with opposite spins are mixed, either time reversal 
symmetry or lattice translation symmetry must be broken 
to form a gapped state with filling fraction 2/m (v = 1/m 
for each spin on average). In Appendix [Fl we prove this 
statement by a careful PSG analysis. 

The analysis in this subsection seems to suggest that 
it is difficult to explicitly construct a FTI wavefunctions 
breaking the S z conservation. In fact, this difficulty is 
due to the formalism of SU(m)^ x SU(m)^ parton con- 
struction. In the next subsection, we propose another 
parton construction formalism which allows us to explic- 
itly write down FTI wavefunctions with spin mixings. 

B. Parton construction of generic TRI FTI states 
in the absence of spin conservation 

In the following we demonstrate that in a new par- 
ton construction formalism, one can write down the elec- 
tron wavefunctions for fully symmetric TRI FTI states 
breaking the S z conservation, and with mean-field terms 
mixing partons with different spins. We introduce the 
following parton construction (m being an odd integer 
and 9 is an arbitrary real number) 

Cf(r) = cos6»n™=i /a,t( r ) + sin6i n™=i3/3,T( r )> 
C4,(r) = -sin6»nT=i/a4( r ) + cos°I\J=i 9P,l( r )- 

(38) 

where / Qj(T and gp i<7 are all fermionic partons each of 
which carries electric charge e/m. It's straightforward 
to see that the electron constructed in this way is in- 
deed a fermion with electric charge e. The TV-electron 
wavefunction at filling fraction v — 2/m (with spin-f 
electrons and N\, = N — N-f spin- 4. electrons) is obtained 
by projection 

AT T JV— iV-f 

$ e ({rl};{rj}) = (0|n^(r!) [] C ^)\ MF ) (39) 
«=i j=i 
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where | M F) is the parton mean-field ground state as will 
be described later. Nf can be any integer between and 
total electron number N. 

The mean-field ansatz can be written as 

(/l CT (r)M Q , Q ,(r l( T|r',a')/ Q '^(r') 
+ 5 t Jr)M' a a , (r, a|r', cr')<V, CT ' (r')) (40) 

where Hamiltonian a /(r, erjr', tr') is the time reversal 
conjugate of M a _ a ' (r, crjr', cr'): 

M; Q ,(r,a|r',a') = aa'Af^r, -cr|r', -a') (41) 

We use spin index cr = ±1 to denote spin f, J,. In the sim- 
plest case when M a ^ (r, cr|r', cr') = (5 QiQ ' Af (r, cr|r', cr'), 
the mean- field ansatz has SU(m) gauge symmetry or 
in other words IGG = SU(m). As discussed ear- 
lier in section [TT1 in a generic case parton mixing 
terms /] ff /^, ff ',a ^ /3 could exist and the IGG of the 
parton mean-field Hamiltonian M ajQ /(r, cr|r', cr') could 
be Z rn , which is the center of group SU(m). Here 
Ma.a'ir, c|r', cr') can be any mean-field ansatz as a so- 
lution of the PSG constraints on the lattice, so that the 
TV-electron wavefunction obtained by projection (|39p pre- 
serves all the lattice symmetries. More precisely, mean- 
field Hamiltonian M QjQ ,/ (r, cr|r', cr') should be invariant 
under a symmetry operation U followed by a SU (m) 
gauge rotation Gry(r, cr) on /-partons. 

For instance, the lattice symmetry groups are shown 
in Appendix [5] for honeycomb lattice and in Appendix 
[C]for checkerboard lattice for spin f / \.. The constraints 
for Gf/(r,cr) will still be those in Appendix [B] on honey- 
comb lattice and Appendix [D] on checkerboard lattice 
in the case when IGG = Z m . Namely Gry(r,t) and 
Gr/(r,!) can be any two solutions of PSG constraints: 
([Blip on honeycomb lattice and (|D13[) on checkerboard 
lattice when IGG = Z m . The symmetry-allowed mean- 
field Hamiltonian M a , a /(r, cr|r', cr') can be obtained in 
the same way shown in Appendix [B] and Appendix [D] 
It's easy to check that spin-mixing terms in mean-field 
Hamiltonian M a , a i(r, a\r' , a') is in general allowed by 
these PSG constraints. 

Meanwhile time reversal symmetry is also preserved 
since the anti-unitary time reversal operation T is real- 
ized by complex conjugation C combined with the follow- 
ing operation 

f a ,*( T ) O cr 3a ._ a (r) (42) 

where spin index cr = ±1 denotes spin t / 4-- One can 
easily check that under time reversal T the spin- 1/2 elec- 
tron operators indeed transform as c a (r) — > crc_ (J (r). Ap- 
parently the above time reversal operation T leaves the 
mean-field ansatz (|40p invariant. At filling v = 2/m (1/m 
filling for each spin on average) we still use the technique 
of inserting flux to enlarge the unit cell by m times. Note 



that when 2ir/m flux is inserted into each unit cell for /- 
partons, an opposite —2ir/m flux must be inserted in each 
unit cell for g-partons to keep the time reversal symmetry. 
Then we fill the lowest m bands for both /-partons and 
g-partons. Each filled band contains N/2 /-partons (or 
g-partons) which correspond to a filling fraction of 1/m. 
It's straightforward to demonstrate that the electron fill- 
ing fraction of state is indeed v = 2/m. There can 
be symmetry-allowed mixing terms between partons with 
different spins in ansatz Af Q Ot /(r, er|r', a'). The mean- 
field state in (|3T)f is a direct product of /-parton state 
and c/-parton state: \MF) = \MF f ) (g) \MF g ). 

There is a real parameter 9 £ [0, 7r/4] which can be con- 
tinuously tuned in our parton construction (|38p . This pa- 
rameter controls the many-body entanglement between 
spin-f and spin- J. electrons in wavefunction (|39[) . It 
should be considered as a variational parameter in vari- 
ational Monte Carlo studies of projected wavefunctions. 
When 9 = clearly there must be equal number of spin-t 
and spin- J. electrons: Nf = N± = N/2 since other com- 
ponents of the many-body wavefunction with N^ =/= N± 
all vanish in ([55)1 . In this case the electron wavefunction 
(|39p is nothing but a direct product of spin-f wavefunc- 
tion $ t(r t) = (0|n^ / 1 2 nLi/ Ct ,t(rl)|A/F/) and spin- 

i wavefunction ^(rj) = (0| ]J^=i Ua=i 9 a ,d^)\MFS). 
This corresponds to the spin-conserved limit when there 
is no entanglement between electrons with different spins. 

When 9 is nonzero, many-body entanglement between 
electrons with different spins is encoded in electron wave- 
function (|39p as long as the spin mixing terms are present 
in the mean-field Hamiltonian (j40|) . And in general the 
component of the many-body wavefunction with an ar- 
bitrary number of spin-f electrons V < < N should 
be nonzero. In a generic case with 9 ^ the many-body 
wavefunction (|39p is complicated and cannot be written 
as a Slater determinant. Now one can see the parton con- 
struction (|38p allows us to write down generic electron 
wavefunctions for TRI FTI states in the absence of spin 
conservation. The spin-conserved TRI FTI wavefunction 
(9 = 0) can be deformed into a generic TRI FTI state 
in the absence of spin conservation {9 ^ 0) by continu- 
ously tuning parameter 9, while keeping the mean-field 
Hamiltonian (|4"U)) unchanged. In the process of tuning 9 
continuously, we expect the low-energy effective theory 
and quasiparticles of such a TRI FTI state to remain the 
same. 

In the end we comment on the low-energy effective 
theory of such a TRI FTI state. In the simplest case 
when M QjQ ,' (r, cr|r', cr') = <5 Q Q /M(r, cr|r', cr'), the mean- 
field ansatz (j4T)|) has a SU(m) x SU(m) gauge sym- 
metry or in other words IGG = SU{m) s x SU(m) 9 . 
Let's assume the filled lowest band of M(r, cr|r', cr') for 
{/a.t/ll a = 1) ' ' ' j TO I partons has a Chern number k. 
Then due to time reversal symmetry the filled lowest 
band for {c/a^/jja = 1, • • • , m} partons will have a Chern 
number —k. And its low-energy effective theory is a 
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SU{m)k x SU(m)_k Chern-Simons theory: 

-■fcepvxTrfaduh + ?bpb v b x ) (43) 

where a M is the ^t/ (m) gauge field coupled to /-partous 
and 6 M is the SU(m) gauge field coupled to g-partons. 
Such a SU(m)f x SU(m) g TRI FTI state will' host non- 
Abclian quasiparticles if k > 1. 

When fc = 1 this is an Abelian TRI FTI state with 
ground state degeneracy m 2 on a torus and anyonic 
quasiparticles. When parton mixing terms /t a fp. a ' , a. ^ 
are present, again the IGG is reduced from SU(m)f x 
SU(m) 9 to Z m * x Z m 9 and the low-energy effective the- 
ory is described by Chern-Simons-Higgs theory, i.e. ef- 
fective action (|4*3"| with a periodic potential due to Bose 
condensation of Higgs fields. Such an Abelian Z„J x Z m 9 
TRI FTI has the same topological properties as an 
Abelian SU(m)l x SU(m) 9 FTI, such as ground state 
degeneracy and quasiparticle charge/statistics. In the 
parton construction ([35]), both SU{m)f x SU{m) 9 and 
Z m * x Z m 9 TRI FTI states are possible candidates for 
a symmetric TRI FTI state in the absence of spin con- 
servation: which state is realized depends on the IGG 
of mean-field amplitudes M aia >(r, er|r', a') and should be 
determined by energetics of wavefunctions (|39p in varia- 
tional studies. 

In the end we comment on the robustness of such 
a symmetry-protected TRI FTI state against perturba- 
tions. As discussed in the previous subsection, in the 
presence of time reversal symmetry and lattice transla- 
tion symmetry, no terms that mix n partons and / 
partons can be added to mean-field ansatz (|40[) . Since 
both /-partons and g-partons form a band insulator, such 
a TRI FTI is stable in the absence of mean-field terms 
mixing /- and g-partons. On the other hand, once time 
reversal symmetry or translation symmetry is broken one 
can always add a mixing term in the mean-field ansatz 
(|40| , which will drive the system from a FTI state into a 
different phase. Therefore the stability of this TRI FTI 
state is protected by time reversal symmetry. 



VI. CONCLUSION 

To summarize, we show that a large class of Abelian 
(and non- Abelian) fractionalized topological phases can 
be constructed on a lattice using parton construction. 
These states preserve all the lattice symmetry and are 
featured by e.g. fractionalized excitations and topologi- 
cal ground state degeneracy. In the spin-polarized case, 
we construct fractional Chern insulator phases belong 
to distinct universality classes even at the same filling 
v = 1/m. Their differences arc characterized by the pro- 
jective symmetry group in the bulk and are protected 
by the lattice symmetry. The low energy gauge groups 
of these states are found to be either SU(m) or Z m . 



Their low-energy physics is described by SU(m)i Chcrn- 
Simons theory and Chern-Simons-Higgs theory respec- 
tively, and they all have m-fold degenerate ground states 
on a torus. We explicitly construct the ground state 
wavefunctions and bulk quasiparticles on the lattice. We 
demonstrate our construction by several explicit exam- 
ples, including non-Abclian FCIs which may be realized 
in a nearly flat band with Chern number C > 1. Fur- 
thermore, we show that when time reversal symmetry 
is present, classes of fractionalized topological insulator 
phases preserving both time reversal symmetry and lat- 
tice symmetries can be constructed. These TRI FTI 
states are characterized by SU (m) x SU (m) or Z m x Z m 
gauge groups. Their electron wavefunctions on the lat- 
tice, which are essentially products of spin-polarized FCI 
states for spin f and its time reversal conjugate, are pro- 
vided. These are stable topological phases even when S z 
conservation is not a symmetry in the electronic system. 
In order to explicitly construct TRI FTI wavefunctions 
with entanglement between opposite spins, we propose a 
new parton construction formalism. It allows one to write 
down generic electron wavefunctions of TRI FTI states, 
which preserve both time reversal and lattice symmetries 
in the absence of spin conservation. Our work provides 
important insight for future numeric study using varia- 
tional Monte Carlo method. 
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After this manuscript was submitted, we note that an- 
other paper— also discusses the technique of inserting 
2 / n/m fluxes in unit cells in order to construct transla- 
tional invariant FCI wavefunctions by parton approach. 

Appendix A: Symmetry group of the honeycomb 
lattice model 

The symmetry group of Haldanc's model is generated 
by the following symmetry operations as shown in FIG. 
[2ja): (1) translations T\ t 2 by Bravais lattice vector Q4.2; 
(2) 7r/3 rotation Cq along z-axis around the honeycomb 
plaquette center; (3) mirror reflection w.r.t. x-o-z plane 
combined with time-reversal operation, labeled as a. 

Note that er is an anti-unitary symmetry since it in- 
cludes time-reversal operation. It acts on the Hamilto- 
nian through a combination of a unitary symmetry op- 
eration and complex conjugation C. Besides, the spatial 
Cq rotation should be accompanied with a corresponding 
spin rotation along S^-axis in a generic electron Hamil- 
tonian with spin-orbit coupling terms. 

We label a lattice site by coordinates (x, y, s) where 
r = xa% + ya2 + r s is its position vector. a\ = a(V3, 0) 
and a,2 = a(\/3, 3)/2 arc two Bravais lattice vectors, s = 
0,1 is the sublattice index with r*o = — a(y/3, l)/2 and 
fi = o(— Vo, l)/2. Under the symmetry operations the 
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(x, y, s) coordinates transform as 

T x : (x,y,s)^(x + l,y,s) (Al) 

T 2 : (x,y,s) -> (x,y + l,s) 

a : (x, y, s) -> (x + y, -y, 1 - s) 

C 6 : (a;, y, s) -> (1 - s - y, x + y + s - 1, 1 - s) 

The multiplication rules of the above symmetry group 
are completely determined by the following algebraic re- 
lations: 

G 6 6 — cr 2 — e, (A2) 
T- l T 2 T{T- x = e, 
T 2 C§T\Cq 1 = e, 
2i 1 CqTiT 2 Gq 1 = e, 

Tf VTio-- 1 = e, 
T^VTjT^V- 1 = e, 
ctCqctCq = e. 

where e represents the identity element of the symmetry 
group. 

Appendix B: Mean-field ansatz of spin-polarized Z m 
FCI states on honeycomb lattice in SU(m) parton 
construction 

1. Projective symmetry group (PSG) analysis of 
different Z m mean-field states 

Since the IGG of a spin-polarized Z m state is 

G e € Z m = { V ■ I mXm \r] = e i2v ^ , k = !,-■■ ,m}(Bl) 

the algebraic relations (|A2[) give us algebraic conditions 
on the gauge transformations Gu(x, y, s): 

G tT {a{v))Gl{r)=r lcr I m , (B2) 
G C6 (C 6 - 1 (r))G C6 (G 6 ~ 2 (r))G C6 (Gf (r)) (B3) 
■Gc 6 (C 6 2 (r))Gc 6 (C 6 (r))G C6 (r) =%/ mxm , (B4) 
G Ti 1 (T 2 - 1 T 1 (r))G T ^(T 1 (r)) 
•G Tl (T 1 (r))G T2 (r) (B5) 
G T ^(T 2 (r))G C6 (T 2 (r)) 
•G Tl (T 1 G 6 - 1 (r))G^ 1 (r) = , C(1 / mxm , (B6) 
G^ i 1 (T 1 (r))Gc 6 (Ti(r))G Tl (G^T^r)) 
■G^ 1 (G 6 - 1 (r))G5 6 1 (r) = r,c 62 / mxm , (B7) 

G 7 ; i 1 (T 1 (r))G CT (T 1 (r)) 
•G^ (T 1 cr- 1 (r))G- 1 (r) = i) ri / mxm , (B8) 
G^(r 2 (r))G CT (r 2 (r))G^(<TT 2 (r)) 

■[G^ 2 1 ( < r(r))]*G- 1 (r) = r^mx™, (B9) 
G CT (r)G^(<x(r)) 
•G;( < rG 6 (r))G C6 (G 6 (r)) (BIO) 



where all the 77's arc Z m quantum numbers {e l2w ~\k = 
1, • • ■ Note that <x is a anti-unitary symmetry so 

it's accompanied with complex conjugate C. We can 
always choose a proper gauge so that Gtj (x, y, s) = 
G T2 (0,y,s) = I mxm . r] Ce i = rjc e 2 = 1 can also be fixed 
by certain gauge choice. After solving the above alge- 
braic conditions we have r/a-i = rja-2 = TI12 and r\ a = 1. 
There arc only two independent Z m quantum numbers 
left, i.e. 7712 and ?7 CT c 6 - In the end we have 

G Tl (x,y,s) =I m xm, G T2 (x,y,s) = ?7f 2 / mxm (Bll) 

G^x,y, s )=^~^ +1)/ V(s), 

Gc 6 (x,y,s) = r)l% 2 9c e (s), 

[5c 6 (0)5c 6 (l)] 3 = m2Vc B , ffo-(l)y*(0) = 1, 
Sc 6 (s).g CT (l - s)y5 6 (s)g*(l - s) = r) aG& - 

where <? CT (s) and gc 6 ( s ) are au SU(m) matrices. 

Different Z m mean-field states are given by gauge- 
inequivalent solutions of equations (|B11|) . 

2. Symmetry conditions on Z m mean- field states 

To choose a representative of all mean-field amplitudes 
M(rjr'), we choose r' = (0,0,0) and label independent 
mean-field bonds as 

[x,y,s] = M{x, y,s|0,0,0) (B12) 

All other symmetry-related mean-field amplitudes can 
be generated by [x, y, s]: e.g. all 1st nearest neighbor 
(n.n.) mean-field amplitudes can be generated by [0, 0, 1] 
through symmetry operations; all 2nd n.n. amplitudes 
can be generated by [0, 1, 0]. 

Using symmetry operations (|A1[) it's straightforward 
to verify the following relations 

T 2 "VG 6 : [-2a;, x, 0] [-2x, x, 0] 

T.f-VGe: [0, x, 0] -> [0, x, 0]+ 
TfVCf : [x, -2x, 0] -> [x, -2x, 0] 
IfVCg: [x,0,0]^[x,0,0]t 
crGfT 1 : [x, x, 0] [x, x, 0] 
Tf r^^crGg- 1 : [x, -x, 0] -> [x, -x, 0]t 

for mean-field amplitudes [x, y, 0] = M(x, y, 0|0, 0, 0) be- 
tween sites 

within the same sublattice and 

T?T 2 - 2x a : [x, -2x, 1] -> [x, -2x, B]t 
crGfT 1 : [x+l,x,l] -> [x+l,x,B] 
r-a»-a 3 w ! +i . C r-a . [_2a;-l,x,l] -> [-2x- l,x,l]+ 
Tf-^Gf: [x,y,l]^[x,y,l]t 

for mean-field amplitudes [x, y, 1] = M(x, y, 1|0, 0, 0) be- 
tween sites from different sublattices. The symmetry con- 
dition (|12|) now becomes constraints on the mean-field 
amplitudes [x,y, s]. 
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In the following we list the symmetry conditions on 
mean-field amplitudes for 1st n.n. u a = [0,0, 1] 

9<r(i)Kgl(Q) = ul, 
.9c 6 (0)ffc 6 (l).9c 6 (0)« a ^ 6 (l)^ 6 (0)^ 6 (l) = ut. 

and 2nd n.n. up = [0, 1,0] 

^(0)^ a (l)tt^(l) 5 t(0) = «J. 



where e represents the identity element of the symmetry 
group. 



Appendix D: Mean- field ansatz of spin-polarized Z m 
FCI states on checkerboard lattice in SU(m) parton 
construction 



Appendix C: Symmetry group of the checkerboard 
lattice model 

The symmetry group of checkerboard lattice model is 
generated by the following symmetry operations as shown 
in FIG.|5[b): (1) translations by Bravais lattice vec- 
tor a'i 2 'i (2) mirror reflection w.r.t. x-o-z plane com- 
bined with time- reversal operation, labeled as P x . (3) 
mirror reflection w.r.t. y-o-i plane combined with time- 
reversal operation, labeled as P y . (4) mirror reflection 
along (x + y)-o-z plane combined with time-reversal op- 
eration, labeled as P xy . 

Again here P x , P y and P xy are all anti- unitary sym- 
metries. They act on the Hamiltonian through a com- 
bination of unitary symmetry operations and complex 
conjugation C. 

Here we also label each lattice site by coordinates 
(x, y, s) where r = xa\ + ya-2 + r s corresponds to its po- 
sition vector, di = a(l, 0) and a,2 = a(0, 1) are two Bra- 
vais lattice vectors and a is the lattice constant. Again 
s = 0, 1 is the sublattice index with r*o = — a(l/2,0) and 
f\ = — a(0, 1/2). Under the symmetry operations the 
[x, y, s) coordinates transform as 

Ti: {x,y,s)^{x + l,y,s) (CI) 
T 2 : (x, y, s) -> (x, y + 1, s) 
Px ■ (x,y,s) ->(x,s- y,s) 

P y - (x, y, s) -> (1 - s - x, y, s) 
Pxy ■ (x,y,s) -> (y,x, 1 - s) 

The multiplication rules of the symmetry group are 
completely determined by the following algebraic rela- 
tions: 

Pl = Pl = Pxy 2 = e, (C2) 

T^TiTxTz 1 = e, 

T 1 - 1 P a T 1 p- 1 = e, 
T^PxT^P- 1 = e, 
T^PyT^Py 1 = e, 

T^PyT.P- 1 = e, 

P\ PxyP2Pxy = 6, 
-^2 PxyP\Pxy 6, 

Px.PyPx.Py 

p p p p — p 
1 xy 1 y 1 xy 1 x — c : 

p p p p — p 
1 xy 1 x 1 xy- 1 y 



1. Projective symmetry group (PSG) analysis of 
different Z m mean-fleld states 

From algebraic relations (|C2[) we have algebraic condi- 
tions on the gauge transformations Gu{x, y, s): 

G^l (x + 1, y - 1, s)G^l (x + 1, y, s) 
■G Tl (x + l,y,s)G T2 (x,y,s) = r) 12 I mxm , (Dl) 
GtI (x, y, s)G Px (x, y, s)G Tl (x, s-y,s) 

■G~pl (X - l,y,s) = TJxlImxm, (D2) 

G T \ {x, y, s)G Pi (x, y, s) [G^ 1 (x,l + s-y, s)] * 

■Gp*(x,y - l,s) = r] x2 I mxm , (D3) 
G^(x,y,s)G Py (x,y,s)[G^(2 -s- x,y,s)]* 

■Gp*(x - l,y,s) = T] yl I mxm , (D4) 
G^l (x, y, s)G Py (x, y, s)G T2 (l-s-x,y, s) 

■Gp*(x,y-l,s) = r)y2l m xm, (D5) 
GtI {x, y, s)G P:cy (x, y, s)G* T . 2 (y, x,l- s) 

■Gpl y (x - l,y, s) = r) xyl I mxm , (D6) 

G^l ( x > y> s ) g p*v ( x > v> s ) G Ti (y> x ' 1 ~ s ) 

■Gp* y (x,y - 1, s) = T] xy 2l m x m , (D7) 
Gp x (x,y, s)Gp x (a;, s-y,s) = rj x I mXm , (D8) 
Gp y (x,y,s)G* Py (l - s - x,y,s) = rj x I mxm , (D9) 
Gp xy {x,y,s)G* P ^(y,x, 1 - s) = r] xy I mxm , (D10) 
Gp T {x, y, s)G* Py (x, s-y, s)G P:r (1 - s - x, s - y, s) 

■G* Py (l - s - x,y, s) = r) Px p y I mXm , (Dll) 
G Px H {x, V, s)G Px (y, x,l- s)G P:ry (y,l-s-x,l- s) 

■ G *P y ( l - s - x iV, s ) = VPxyPxtmxm- (D12) 

where again all 77's are Z m quantum numbers {e l27r ~ \k = 
1, • • • ,m}. Note that P x , P y , P xy are all a anti-unitary 
symmetries and they are accompanied with complex con- 
jugate C. We can always choose a proper gauge so that 
G T2 ( x ,y, s ) = G Tl (x,0,s) = I mXm - Vxyi = 1 can also 
be fixed by certain gauge choice. After solving the above 
algebraic conditions one can sec rj y i = r) X 2 = Vx = Vy ~ 

Vxy = 1, Vxy2 = Vxyl = 1, Vy2 = Vxl ■ lt turnS 0ut 

V12, Wxi, VPxPyi VP xy Pa: arc tne on ly ^ 0UT independent 
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Z m quantum numbers. In the end we have 

Gt 2 (x, y, s) = I mxm , G Tl (x, y, s) = r]f 2 I mx jpi3) 

Gp v (x, y, s) = Vxi9y(s), G Pxy (x, y, s) = rf^g xy {s), 
9x(s)g*(s) = 9y{s)g* y {s) = 9xy(s)g* y {l - s) = 1, 

[9x(s)g*(s)] =Vxi r )P*P y , 
9xy{s)g* y {l - s)g xy (l - s)g*(s) = ?y PxH P x . 

where g x (s), g y (s) and g X y{s) are all SU(m) matrices. 

Different Z m mean-field states are given by gauge- 
inequivalent solutions of equations (|D13[) . 

2. Symmetry conditions on Z m mean- field states 

A representative of all mean-field amplitudes M(r\r') 
is given by r' = (0, 0, 0) and we label independent mean- 
field amplitudes as 

[x,y,s] = M(x, y,s|0,0,0) (D14) 

All other symmetry-related mean-field amplitudes can 
be generated by [x,y,s}\ e.g. all 1st nearest neighbor 
(n.n.) mean-field amplitudes can be generated by [0, 0, 1] 
through symmetry operations; all 2nd n.n. amplitudes 
can be generated by [0, 1,0] and [1,0,0] etc. . 

Using symmetry operations (|C1|) it's straightforward 
to verify the following relations 

Tl~ x TlP x P y : [x,y,0]-+ [x,y,0] t 
IfP x : [0,x,0] ->■ [0,1, 0]t 
T^Py: [0, x, 0] -> [0, x, 0] 
Tf^Pyi [x,0,0] -> [x,0,0]t 
P x : [a, 0,0] -> [a, 0,0] 

for mean-field amplitudes [x, y,0] = M(x, y, 0|0, 0, 0) be- 
tween sites within the same sublattice and 

TfT^ x P xy : [x, -x, 1] -»■ [x, -x, l]t 
TfT£P x P y P xy : [x, x + 1, 1] -> [x, x + 1, l]t 

for mean-field amplitudes [x, y, 1] = M(x, y, 1|0, 0, 0) be- 
tween sites from different sublattices. The symmetry con- 
dition (|T2|) now gives us constraints on the mean-field 
amplitudes [x,y, s\. 

In the following we list the symmetry conditions on 
mean-field amplitudes up to 3rd nearest neighbor. For 
1st n.n. u a = [0, 0, 1] 

9P,M<9k v W=<- 
for 2nd n.n. up x = [1,0,0] 

■q x ig x {Q)u* px gl(Q) = up x , 
9y(0)u* Px gl(0) = u\ x . 



for 2nd n.n. up y = [0, 1, 0] 

VxiVi2 9y(Q)u} v 9l(0) = u Py , 
g x (0)u* py gt(0) = u\ y . 

for 3rd n.n. w 7 = [1, 1,0] 

4=^(0)^(0)^(0)4(0) (D15) 

Appendix E: Ground state degeneracy of SU(m)\ 
Chern-Simons theory on a torus 

In this section we calculate the ground state degener- 
acy of a SU{m)x Chern-Simons theory 

£su(m)i = ^nvX^{i'ii9 v ax + ^a^ax), e SU(m) 

(El) 

on a torus with m being an odd integer. In Ref. |45| the 
ground state degeneracy of SU(3) q parton states on a 
torus is calculated. We follow their strategy^ and gen- 
eralize it to the case of SU{m)\. 

After choosing a gauge ao = 0, the classical config- 
uration of gauge fields a M are constrained by the fol- 
lowing condition b = e^didj = 0. The gauge field 
configuration is fully determined by Wilson loop opera- 
tors Ui.2 = V\e ^ C1 ' 2 ""^"l for the two non-contractible 
loops ci t 2 along x\ and Xi directions on the torus. 
These two loop operators commute with each other since 
U\UiV\u\ — 1 is the Wilson loop operator for a con- 
tracture loop. Therefore by a global SU{m) transforma- 
tion, we can see lie in the maximal Abclian subgroup 
of SU(m), which is generated by the Cartan subalgcbra 
of Lie algebra su{m). We choose the generator of its 
Cartan subalgebra to be m x m matrices 

ff 1 =Diag(l > -l > 0,--. > 0), (E2) 
g 2 =Diag(l,l, -2, (),■■■ ,0), 
g 3 = Diag(l, 1,1,-3,0,- •• ,0), 

g™- 1 =Diag(l,-.- ,1,-m-l). 
Therefore the gauge field configuration is given by 

m—l 

a Q = 0, aj = ^2 a j9' '■ 
1=1 

Plugging this into (|E1[) we have 

1=1 

In this way the original Chern-Simons theory of non- 
Abelian SU(m) gauge fields is reduced to (m — 1) dif- 
ferent Chern-Simons theory of U(l) gauge fields a^, I = 
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1, ••■ ,m — 1. In addition to these Abelian U(l) gauge 
structures, there are discrete SU(m) gauge transforma- 
tions generated by Wi € SU(m) which leaves this max- 
imial Abelian subgroup invariant. The low-energy de- 
grees of freedom are describe d 51 ' 52 by vectors u(t) and 
v(t) with 



2tt 



2n 



a{(xi,x 2 ,t) = — a I 2 (x 1 ,X2,t) = —vi(t). (E3) 



Li 



L, 



on a L\ x L 2 torus. The effective action becomes 



C 



m — 1 

E 

i=i 



2tt > '1(1 + l)wjitj 



(E4) 



This immediately leads to the following canonical com- 
mutation relations 



[ui,Vj 



1 



2tt 1(1 + 1 



(E5) 



In other words, the conjugate momenta of coordinate it/ 
is 27r/(/+l)u/. Due to uncertainty principle, they cannot 
be fixed simultaneous in a quantum state. 

There are large gauge transformations {£// = 
exp (Yl-KXig 1 jLi}, i = 1,2} acting on the partons, which 
leaves both the electron operators and the physical elec- 
tron states invariant. The conjugate variables transform 
in the following way under the large gauge transforma- 
tions 

ui ~ ui + 1, vi ~ Vi + 1. 

therefore they live on a torus of size lxl. As a result 
the wavefunctions can be written as 



4>{u) = ) 



E c " p27 



(E6) 



where n is a (m— l)-dimcnsional vector of integers. Since 
the conjugate momentum of ui is 2irI(I + l)vi we have 



ijj(u) 



E 



7=1 



(E7) 



The large gauge transformation uj ~ Uj + 1 enforces 
the periodic condition en = cn+sn where Sftj = 1(1 + 
l)Zi, Zi G Z. 

Aside from the large gauge transformations U* , there 
are other large gauge transformations under which 



/ ui 

\Um-l/ 



(E8) 



-1 


1 


1 •• 


1 





-2 


1 • • 


1 








-3 •• 


1 








•• 


■ —m — 1 




and similar relations hold for v. Here {di, I = 2, ■ ■ • , m} 
are arbitrary integers and d\ = — ^2i=2 The partons 
transform as f a (xi 1 X2) ~^ exp(i2nxid a / Li)f a (xi, x 2 )- 
As mentioned earlier, in addition to the large gauge 
transformations, there are discrete SU(m) gauge trans- 
formations essentially generated by permutation between 
the m species of partons. For example, the permutation 
Pi 2 transforms the partons as f\ O f 2 and thus we have 



Pis 



u = (ui,u 2 , 
v= (v 1} v 2 , ■ 



(~Vl,V 2 ,- ■ -) T 



There are in total m(m— 1)/2 permutations {P Q( g|a < /?} 
that leaves the maximal Abelian subgroup of SU (m) in- 
variant. All these gauge transformations give us further 
equivalence conditions on u, v, and thus on coefficients 
Cft. In the end, the number of independent coefficients Cft 
equals the ground state degeneracy of SU(m)i Chern- 
Simons theory on a torus. With the help of computer 
programming it's easy to check the ground state degener- 
acy of SU(m)i Chern-Simons theory on a torus is m, and 
they correspond to the following (gauge-incquivalent) 
gauge field configurations 



vi = 0, (J = l,2,--- ,m-2); 

fm-i = i fc=l,2,---,m. (E9) 



Physically they correspond to inserting 2irk/m flux 
through the hole along x\ direction. An explicit example 
of SU(3)i Chern-Simons theory is given in Ref. liH 



Appendix F: Symmetry analysis of 
time-reversal-invariant FTI states 



In this section we use projected symmetry group (PSG) 
to analyze symmetric mean-field ansatz realized in a 
time-reversal invariant lattice model in SU (m) parton 
construction. We require the corresponding electron 
state to preserve both time reversal symmetry T and all 
lattice symmetries. 

In the following we use the time-reversal-invariant 
checkerboard lattice mode l 11 ' 16 as an example. The dis- 
cussion is completely general and can be applied to any 
lattice model. The symmetry group is generated by lat- 
tice symmetry (jCl[) together with time reversal symme- 
try: 



(x,y,s,a) (x,y,s,l - a) 



(Fl) 



where a = 0/1 represent the electron spin t / 4- Note 
that here P x , P y and P xy are the combination of spatial 
7r rotations (along x, y and x + y-axis) and spin rota- 
tions. So all these three generators flip spin: a — > 1 — a. 
They are not accompanied with time-reversal operation 
anymore and they are now unitary symmetries. Time re- 
versal itself is still an anti-unitary symmetry realized by 
complex conjugation operation C. The algebra structure 
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of this symmetry group is determined by (IC2[) together 
with 

T 2 = e; (F2) 
g-'T-'gT = e, 

.9 = T1.2, P x , Py, Pxy 

Consider a IGG = Zj n x mean-field ansatz and 
we have the following algebraic conditions for the 
gauge transformations Gu(x, y, s, a) associated with U = 
T,T lj2 

G^(x,y,s,<7)G T (x,y,s,a)G^ i (x,y,s, I- a) 

G^, 1 {x-l,y > s > a)=riTi, (F3) 

GtI ( x > V> s > &)Gt{x, y, s, a)G* T2 (x, y,s,l- a) 

G T 1 (x,y-l,s,<j)=? lT2 . (F4) 

since T-f^TTi^T -1 = e. Here tjti, VT2 are all Z m quan- 
tum numbers {e i2j ^\a = 1, • • • , m}. As shown in Ap- 
pendix [HI one can always choose a certain gauge so that 

GT 2 {x,y,s,a) = I m xm, Gt^ (x, y, s, a) = r]i 2 {cr)I m x m - 

where 7712(f), c = 0/1 are also Z3 quantum numbers. 
From (|F3p and (|F4j) one can solve out 

G T (x,y,s,a) = Vti11t29t{s,o), 

vMvM = 1. (F5) 

Therefore when IGG = Z^ x Z^ for the parton mean- 
field state, we can simply choose 7712 (4-) = V^it) to pre- 
serve the time-reversal symmetry. However, when there 
are spin mixing terms in the mean-field parton Hamil- 
tonian, the gauge symmetry is broken further down to 
IGG = Z m and we must have 7712(f) = 7712(f). From 
(|F5[) we know in order to preserve time-reversal sym- 
metry, we must have 7712(f) = 7712(f) = 1 m the case 
IGG = Z m . As a result there is no way to insert 27r/m 
flux into each unit cell in the parton mean-field ansatz 
if we require time-reversal symmetry and lattice trans- 
lation symmetry. Therefore the technique of enlarging 
unit cell by m times in order to construct a gapped state 
with filling fraction v = 1/m failed. This suggests that 
a gapped FTI state with gauge structure IGG = Z m at 
filling fraction 1/m would break cither time reversal or 
lattice symmetries. 

Appendix G: The honeycomb lattice model: four 
different examples of spin-polarized FCI states at 

v = 1/3 

Haldane's model on honeycomb lattice^ has been 
shown to support nearly flat bands with non-Zero Chern 
numbers^. Its symmetry group is shown in Appendix fAl 
and FIG. [2J We label a lattice site by coordinate (x, y, s) 
as explained in Appendix [A] By SU{3) parton construc- 
tion, we use PSG to classify different spin-polarized Z3 




FIG. 9: (color online) Mean-field ansatz of the parent SU(3) 
parton states associated with FCI states HC1 and HC3 on 
honeycomb lattice. The solid line for NN bonds represents 
real hopping amplitude a. NNN bonds represents complex 
hopping amplitude (3 along direction of the arrow. The triple 
arrow means the original hopping amplitude along its direc- 
tion should be multiplied by a phase factor 7712. Note in the 
mean-field ansatz the lattice translation along 0,2 direction is 
explicitly broken by flux insertion and the unit cell is tripled. 
However the lattice translation symmetry is preserved in the 
corresponding electron states after projection (JSJl. 

mean-field ansatz as shown in Appendix [B] In the fol- 
lowing we show two mean-field states belonging to differ- 
ent universality classes. They correspond to two gauge- 
incquivalent solutions to (|F31 1[) . Their two SU(3) parent 
states also have different PSGs and thus distinct. We 
show mean-field amplitudes up to NNN in FIG. M 

1. FCI state HC1 with a xy = 1/3: the Z 3 state and 
its parent SU(3) state 

In spin-polarized Z 3 FCI state HC1 the gauge trans- 
formations Gu{x,y, s) associated with lattice symmetry 
U are listed below: 

G Tl ( x >y> s ) = hx3, G T2 (x,y,s) = 77^/3x3, 
<^o-{x,y, s) — ?7 12 ^3x3, 

G, \ xy+x(x — l)/2 T 

c 6 {x,y,s) = ?7 12 7 3>< 3 

where 7712 = cxp(i27r/3). As shown in Appendix 151 the 
symmetry allowed mean-field amplitudes are: 

(I) For nearest neighbor (NN) amplitude u a = 
M(0,0, 1|0, 0, 0) 

u a = = it* (Gl) 

in other words u a can be any real symmetric 3x3 matrix. 
All other NN amplitudes can be generated from u a by 
symmetry operations through (|12[) . 
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HC1: parent SU(3) state 
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HC2: Z3 state 
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FIG. 10: (color online) The parton band structure of the 
parent SU(3) ansatz (|25p of HC1 state. Each band is 3- 
fold degenerate, corresponding to 3 parton flavors f 1,2,3- We 
plot the dispersion along 1 -> A' -> A/ -> 1 as shown in 
FIG. |3{a). Hopping parameters are chosen as a = — 1 for 
NN, = -0.22exp[0.2i] for NNN. The Chern numbers of the 
6 bands are {1, —2, —2, 1, 1, 1} in a bottom-up order. 



FIG. 11: (color online) The parton band structure of a Z3 
mean-field ansatz: HC2 state. There are 18 non-degenerate 
parton bands. We plot the dispersion along T — > K — > M — > T 
as shown in FIG. EJa). The Chern numbers of the lowest 6 
bands are {1, 1, 1, —2, —2, —2} in a bottom-up order. 



below: 



(II) For next nearest neighbor (NNN) amplitude up 
M(0, 1,0|0,0,0) 



up = uj 



(G2) 



i.e. up can be any symmetric 3x3 matrix. All other 
NNN mean-field amplitudes can be generated from u a 
by symmetry operations through (|12p . 

Its parent 517(3) mean- field state includes both NN 
and NNN hopping terms u a — 0/3x3, up = /3i3x3, as 
shown in FIG. |H1 With hopping parameters chosen as 
a = — 1, /3 = — 0.22 cxp[0.2i], the parton band structure 
is shown in FIG. [TDJ In the bottom-up order the Chern 
numbers are {1,-2,-2,1,1,1} for the six bands (each is 
3-fold degenerate) and these band structures persist in 
a large parameter range. As we add small terms mixing 
different species of partons to break the gauge symmetry 
down to Z 3 , we have 18 non-degenerate parton bands and 
the Chern number for the lowest three bands are all +1. 
And the lowest three parton bands are well separated 
from the other 15 bands by a large gap. As a result 
the Z 3 FCI state we obtained by filling the lowest three 
bands has a Hall conductance 3 • (1/3) 2 = 1/3 in the unit 
of e 2 /h since each parton carries electric charge e/3. 



2. FCI state HC2 with o xy = 1/3: the Z z state and 
its parent 517(3) state 

In Z 3 FCI state HC2 the gauge transformations 
G[/(x, y, s) associated with lattice symmetry 17 are listed 



G Tl (x,y,s) = I 3 x3, G T2 (x,y,s) = rjf 2 I 3x3 , 
^a\X, y, s) — r] 12 + 3x3, 

Gc 6 (a;, y, 0) =7?^ v ;/ / 3X 3, 



G Ce (x,y, 1) = V12 



xy-\-x(x — l) jl 



Diag{ 



,i27r/9 „i27r/9 



-i47r/9 



where 7/12 = cxp(i27r/3). 

As shown in Appendix [B] the symmetry allowed mean- 
field amplitudes are: 

(I) For nearest neighbor (NN) amplitude u a = 
M(0,0, 1|0, 0, 0) 



e i7r / 9 au 



e -iir/18 



iir/9 



Q-22 
i7r/18. 



— i7r/18~. 

e ' ai3 

i7r/18~, 
e ' Qf23 



a 13 e l "'^a 23 e i27r/9 a 3 3. 



(G3) 



where all ctij are real parameters. All other NN ampli- 
tudes can be generated from up by symmetry operations 
through (|T2"j), 

(II) For next nearest neighbor (NNN) amplitude up = 
M(0,1,0|0,0,0) 



up 



P11 
/81a 



-in/3 



Pl2 
P22 



323 e 



-iir/3 



/3i 3 e i7r/3 ' 
/? 23 e i7r/3 



J 3-'S 



(G4) 



where all are complex parameters. All other NNN 
mean-field amplitudes can be generated from up by sym- 
metry operations through (|12[) . 

Its parent 517(3) mean-field state include only NNN 
hopping terms, i.e. u a = and up = /3I 3x3 . It is demon- 
strated by the NNN bonds shown in FIG. Appar- 
ently this mean- field ansatz is very different from HC1 
state since its parent 517(3) state doesn't allow NN hop- 
ping terms. We choose the following hopping parameters: 
an = a 2 2 = 0:33 = 2, a?i 2 = 0.07, ai 3 = 0.04, a 2 3 = 
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FIG. 12: (color online) The parton band structure of HC3 
state, a non-Abelian 517(3) FCI state on honeycomb lattice. 
Each band is 3-fold degenerate, corresponding to 3 parton 
flavors /i,2,3- We plot the dispersion along F — s- K — > M — > V 
as shown in FIG. [3ja). Its mean-field ansatz are shown in 
FIG. [9] with 7712 = exp(— i27r/3). Hopping parameters are 
chosen as a = -1 for NN and /3 = -0.22 exp[0.2i] for NNN. 
The Chern numbers of the 6 bands are {2, 2, —1, —4, 2, —1} in 
a bottom-up order. 
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FIG. 13: (color online) The parton band structure of CB3 
state, a non-Abelian SU (3) FCI state on checkerboard lattice. 
Each band is 3-fold degenerate, corresponding to 3 parton 
flavors /i,2.3- We plot the dispersion along T —¥ K —¥ M —¥ 
r as shown in FIG. [3jb). Its mean-field ansatz are shown 
in FIG. with 7712 = exp(i27r/3). Hopping parameters are 
chosen as a = e in/12 for NN, f3 x = 0.6 = j3 y for NNN and 
7 = 0.5 for NNNN. The Chern numbers of the 6 bands are 
{2, 2, —1, —4, 2, —1} in a bottom-up order. 



0.05 and /3 n = hi = A33 = 0.1 exp[0.2i], /3 12 = 
0.03exp[i0.7], /3 13 = 0.07 exp[i0.2], /3 23 = 0.06 exp[i 1.9]. 
We find that in a very large parameter range the Chern 
number of the lowest three parton bands remains to be 
+1 and they are separated from the other bands by a 
large gap, as shown in FIG. [TT] By filling these 3 lowest 
bands we obtain the Z 3 FCI state HC2 which also has 
Hall conductivity 1/3 in the unit of e 2 /h. 



Appendix H: Two examples of non-Abelian 
spin-polarized ST/ (3) FCI states 

As discussed in section ITVl by partially filling a topo- 
logical flat band with Chern number C > 1, it might be 
possible to realize non-Abelian FCI states. Here we show 
two examples of spin-polarized SU(3) FCI states with 
v = 1/3, one (labeled as HC3) on honeycomb lattice and 
the other (labeled as CB3) on checkerboard lattice. In 
particular these states are likely to be non-Abelian FCI 
states associated with SU(3)2 Chern-Simons theory^ 
since the each parton species fills a lowest band with 
Chern number +2. Their Hall conductance is <r xy = | 
in the unit of e 2 /h. By braiding non-Abelian quasiparti- 
cles in these FCIs one can carry out universal quantum 
computations^^. 



1. Honeycomb lattice: 577(3) FCI state HC3 with 

a xy — 2/3 and v = 1/3 

The SU(3) FCI state HC3 is an example of non- 
Abelian states on the honeycomb lattice. It includes both 
NN and NNN hopping terms u a = 0/3x3 , up = ftl 3 x 3 



in its mean-field ansatz, as shown in FIG. with 7/12 = 
cxp(— i27r/3). It is different from the parent SU(3) state 
of HC1, since in HC1 we insert — 27r/3 flux per unit cell 
and here we insert 27r/3 flux in HC3 state to enlarge the 
unit cell. With hopping parameters chosen as a = 1, 
f3 = 0.35 exp[0.2i], its parton band structure is shown 
in FIG. Q21 In the bottom- up order the Chern num- 
bers are {2,2,-1,-4,2,-1} for the six bands (each is 
3-fold degenerate) and these band structures persist in a 
large parameter range. And the lowest parton band is 
separated from the other 5 bands by a sizable gap. As 
a result the SU (3) FCI state HC3 we obtained by filling 
the lowest 3-fold-degenerate band has a Hall conductance 
&xy = 3 ■ 2 ■ (1/3) 2 = 2/3 in the unit of e 2 /h, since each 
parton carries electric charge e/3. 



2. Checkerboard lattice: SU(3) FCI state CB3 with 

a xy = 2/3 and v — 1/3 

The SU(3) FCI state CB3 is an example of non- 
Abelian FCI states on checkerboard lattice. It in- 
cludes NN hopping u a = a/3 X 3, NNN hopping u ax = 
e- i7T / 3 f3 x I 3x3 , u ay = /3 y I 3x3 and NNNN hopping u 7 = 
7/3x3 in its mean-field ansatz, as shown in FIG. [7] with 
i]i2 = cxp(i27r/3). It is different from the parent SU(3) 
state of CB2, since in CB2 27r/3 flux are inserted into 
each unit cell while here we insert — 27r/3 flux for CB3. 
Choosing parameters a = e 17r / 12 , (3 X = 0.6 = j3 y and 
7 = 0.5 we have the Chern numbers of the 6 bands for 
each parton species as {2,2,— 1,-4, 2,-1} and the low- 
est band is well separated from other bands, as shown 
in FIG. [T3] This qualitative band structure persists for 
a large parameter range. Each band of a SU (3) parton 
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ansatz is 3-fold degenerate, corresponding to the 3 par- state CB3 whose Hall conductivity is a xy = 2/3 in the 
ton flavors /i 2,3- By filling the 3-fold-degenerate lowest unit of e 2 /h. 
band (with Chcrn number +2) we obtain a SU(3) FCI 
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